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I.  INTRODUCTION 


In  the  last  year  the  national  program  of  particle  beam  drivers  for 
Inertial  Confinement  Fusion  (ICF)  has  been  shifted  from  Relativistic  Electron 
Beams (REB)  to  light  ion  beams.  The  major  reasons  for  preferring  Ion 
beams  Mere  the  high  efficiency  of  generating  them  using  present  pulsed  power 
generators  coupled  to  Ion  diodes,  focusing  and  transporting  the  beams,  and 
bunching  them  In  plasma  channels.  Many  of  these  results  have  been  achieved 
by  the  NRL  program  supported  by  JAYCOR.  In  particular.  Me  have 
proposed  the  possibility  of  Increasing  the  power  of  the  Ion  beam  by  reducing 
the  pulse  length  from  the  100  ns  to  20  or  10  ns  needed  for  ICF  targets. 


In  the  present  report  concentrate^  the  physics  of  propagating 
Intense  Ion  beams  In  plasma  channels.  In  particular,  we  have  studied  the 


theoretical  aspectsAof  the  possible  effects  during  Ion  beam  transport  on 
the  structure  of  electromagnetic  fields  In  the  plasma.  The  Ion  orbits 
were  solved  for  different  classes  of  plasma  channels  Including  bumpy  ones 
in  order  to  evaluate  the  sensitivity  of  beam  losses  due  to  channels  that 
deviate  from  Idealistic  shapes.  The  consideration  of  different  plasma 
instabilities,  the  hydrodynamic  motion  of  the  plasma,  propagation  In  fila¬ 
mentary  channels  and  beam  energy  losses  put  realistic  limits  on  how  much 
Ion  current  can  be  propagated  using  the  simple  steady  plasma  channel  concept. 
Axial  bunching  of  the  beam  (via  time  dependent  voltage  ramping  of  the  ion 
diode)  was  analyzed  and  found  to  be  a  sound  concept  as  long  as  no  more  than 
an  order  of  magnitude  pulse  compression  Is  attempted.  The  radial  compression 
Is  viable  on  a  small  scale  In  tapered  channels  (giving  a  factor  of  2  Increase 
In  current  density)  and  on  a  large  scale  if  self  magnetic  fields  fully  pene¬ 
trate  the  media  (giving  rise  to  an  order  of  magnitude  Increase  In  Ion  beam 
current  density).  Further  theoretical  Investigation  on  the  details  of  the 


bove  mechanisms  are  presently  In  progress  including  new  plasma  channel 
schemes  that  help  to  push  upwards  the  limits  on  the  ion  currents. 
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II.  SUMMARY  OF  WORK 


The  theoretical  studies  were  performed  along  the  same  approach  as 
envisioned  In  the  proposal  with  additions  that  further  clarify  each  subject. 

1.  In  the  area  of  beam-plasma  Instability,  the  following  results  were 
obtained: 

a.  The  Initial  plasma  current  was  found  to  be  about  50  kA  for 
present  Ion  beams.  The  gas  deuterium  was  shown  to  be  preferential 
to  hydrogen  with  the  needed  density  10”5  g/cm3. 

b.  The  electromagnetic  Welbel  and  Whistler  modes  were  studied 
and  found  to  be  unstable;  however,  they  do  not  grow  fast  enough 
to  seriously  affect  beam  transport. 

c.  In  the  transport  channel,  where  electrostatic  modes  are 
colllslonally  stabilized,  these  electromagnetic  modes  are  the 
most  dangerous  Instabilities.  But  for  lower  plasma  densities 
and  higher  temperatures  the  electrostatic  inodes  are  the  fastest 
growing. 

2.  Ion  Beam  Bunching. 

a.  The  Ion  beam  energy  fluctuations  cause  pulse  broadening  but 
power  multiplication  by  a  factor  of  5  is  easily  obtainable. 

b.  Ion  orbits  in  filamented  channels  have  been  studied  resulting  in 
limits  on  the  tolerable  level  of  jitter  in  the  Ion  beam  current 
density  profile.  If  these  limits  are  violated,  beam  transport 

will  be  affected,  particularly  at  the  tall  of  the  beam. 
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c.  Exploring  the  phase  space  picture  of  the  ion  orbits  as  defined 
at  the  initial  injection  point  leads  to  conditions  for  minimizing 
the  amount  of  beam  rippling  in  the  channel.  Voltage  ramping, 

ion  scattering  during  transport  and  non-ideal  magnetic  field  pro¬ 
files  all  decrease  the  usual  ripple  phase-mixing  distance. 

d.  The  optimum  plasma  and  beam  parameters  for  transport  of  a  given 
current  density  ion  beam  have  been  identified  to  have  injection 
angles  a  <  0.2  (values  of  a  *  0.1  are  preferable  but  difficult  to 
achieve)  and  ion  currents  of  less  than  1  MA  in  the  1  cm2  channels 
studied.  The  major  limit  on  bunching  is  not  abut  voltage  fluctua¬ 
tions. 

3.  Radial  Compression  of  Ion  Beams. 

a.  Studying  ion  orbits  in  a  tapered  z-discharge  plasma  channel 
has  indicated  that  radial  compression  of  the  beam  used  in  con¬ 
junction  with  bunching  can  enhance  bunching  power  multiplication 
by  at  most  a  factor  of  1.3. 

b.  The  ion  orbits  computed  in  the  steep  gradient  magnetic  fields 
(generated  by  either  imploding  plasma  channels  or  expanding  ones) 
are  similar  to  those  of  constant  current  profiles  in  the  plasma 
but  need  somewhat  larger  current  for  confinement.  The  major 
result  is  that  one  has  to  start  with  initial  current  above  the 
initial  minimum  current  required  for  confinement,  since  plasma 
expansion  sweeps  out  the  linear  magnetic  field  with  it  generating 
steep  gradients. 

c.  The  self  fields  and  Ion  orbits  were  determined  for  the  case 
where  current  neutralization  does  not  occur.  Focusing  at  a 
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distance  shorter  than  a  quarter  of  a  betatron  wave  length  was 
found.  The  self  focusing  increased  the  current  density  by  more 
than  a  factor  of  ten. 

d.  The  sheet  Ion  beam  converging  in  cylindrical  geometry  with  a 
center  line  of  symmetry  and  a  plane  perpendicular  to  the  line  with 
anti -synmetri c  magnetic  fields  was  shown  to  follow  a  similar  Ion 
orbit  pattern  to  the  tapered  z-discharge  (3a)  but  with  an  additional 
current  density  Increase  of  r‘^.  It  nay  prove  useful  in  final 
focusing  near  the  pellet. 


! 
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III.  THEORY  OF  ION  BEAM  BEHAVIOR  IN  PLASMA  CHANNELS 

We  first  discuss  the  beam-plasma  Instabilities  in  a  given  plasma 
channel  describing  our  basic  considerations  for  choosing  the  specific  set 
of  Initial  plasma  parameters. 

1.  Beam-Plasma  Instabilities. 

The  Initial  values  of  the  MHO  fields  are  dictated  by  the  fact  that 

the  Ion  beam  was  assumed  to  have  about  1  cm2  cross  section  that  is  a  little 

smaller  than  the  confining  low  density  plasma.  The  plasma  was  chosen  to  be 

hydrogen  (or  deuterium  that  has  twice  the  mass  density  for  the  same  plasma 

density  number  of  electrons)  and  Is  fully  Ionized  In  the  channel.  The  needed 

plasma  current  Is  about  50  kA,  enough  to  confine  5  MeV  protons  with  spot 

size  /Z  smaller  than  the  plasma  channel  and  injection  angle  of  0.15  radians. 

Since  the  plasma  stopping  power  Is  (when  fully  ionized)  lO3/^  MeV/g/cm2 , 

where  V^lon  energy  In  MeV,  and  since  propagation  over  distances  of  5  meters 

18  -3 

may  be  necessary  hydrogen  densities  less  than  4  x  10  cm  should  be  used 

for  less  than  10%  energy  loss.  At  the  same  time  the  hydrogen  density  cannot 

17  -3 

be  decreased  by  large  numbers  to  say  10  cm  since  the  plasma  expansion  is 
too  rapid  for  ion  beam  current  densities  of  *  0.5  MA/cm2.  The  plasma 
expansion  Is  computed  from  Jp  x  8,  the  pressure  induced  by  the  return  current 
Op*  -  Jj.  This  gives  an  equation  for  the  radial  velocity  Vr, 

dV„  x  B  x  6 

r  _  d  1 

p  Ht - c  -  "1 —  ’ 

where  p  Is  the  plasma  density  (g/cm3).  For  a  linear  estimate 
0,  Bt 

1/  _  _ J _ 
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which  Implies  that  a  1017  cm’  hydrogen  plasma  reaches  a  velocity  of  10* 

cm/s  during  t  *  50  ns  for  B  *  2  x  10"  Gauss  (I  *  50  kA,  r  =  0.5  cm).  This 

expands  the  plasma  radius  beyond  2  cm  and  thus  ion  beam  current  density 

reduction  by  a  factor  16  will  occur.  In  addition,  very  large  return  current 

electric  fields  due  to  plasma  motion  (if  a  -*■  ®  the  electric  field  is  given 
Vp  x  B 

by  E  *  -  — - —  >  10"  V/cm)  cause  unacceptable  energy  losses  (all  5  MeV  are 
lost  in  5  meters  of  transport  for  the  above  example).  Thus,  hydrogen 
density  must  exceed  1  x  1018  cm"3. 

An  initial  plasma  temperature  of  a  few  eV  is  obtained  in  equilibrium 
for  a  50  kA  channel  current.  The  heating  of  the  channel  electron  to  50  eV 
is  due  to  the  plasma  stopping  power  (resistive  heating  is  negligible).  The 
ion  beams  were  thus  studied  in  the  above  parameter  regimes  taking  into 
account  the  effect  of  different  temperatures  and  different  current  profiles 
as  the  plasma  undergoes  expansion.  The  plasma  responds  to  the  ion  beam  on 
a  long  hydrodynamic  time  scale1  and  on  a  short,  velocity  phase-space  time 
scale  that  can  excite  fast  field  instabilities.  We  assume  that  the  two  time 
scales  are  separated  enough  that  we  may  treat  them  with  a  semi -static  approach, 
computing  a  Vlasov  formalism  during  different  hydrodynamic  states.  In 
future  work  we  may  resolve  the  coupling  between  the  two  phenomena  when  no 
large  time  separation  exists. 

Velocity  space  electromagnetic  instabilities,  which  lead  to  azimuthal 

and  radial  current  bunching  in  the  channel,  have  been  studied.  In  particular, 

growth  rates  for  the  Weibel  (k.*B  *  0,  k.*V2»0)  and  Whistler  (k  x  B  *0, 

k;^  *  0)  Instabilities  are  derived  In  Appendix  1.  Although  electrostatic 

Instabilities  typically  have  larger  growth  rates  than  electromagnetic 

Instabilities,  here  the  Weibel  and  Whistler  modes  have  the  fastest  growth 

2 

because  of  colllslonal  damping  of  the  electrostatic  modes. 
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The  plasma  electron  current  is  the  main  driver  for  the  Weibel  insta¬ 
bility  and  the  ion  beam  current  is  the  main  driver  for  the  Whistler  insta¬ 
bility.  Both  instabilities  have  relatively  slow  growth  rates  so  that,  at 
most,  1-2  e-folds  occur  during  beam  passage  for  typical  systems.  Nonlinear 
saturation  probably  is  not  reached  during  the  beam  propagation  time,  however, 
the  late  time  level  of  current  bunching  can  be  estimated  using  the  predicted 
e-folding  numbers  and  the  Initial  perturbations  in  the  system.  Such  pertur¬ 
bations  are  most  likely  produced  by  jitter  in  the  beam  current  density  (see 
Appendix  2).  Since  both  instabilities  are  nonconvecting,  beam  ions 
generated  at  the  tail  of  the  pulse  are  most  affected  by  the  current  bunching. 

2.  Ion  Beam  Bunching. 

The  effect  of  diode  voltage  fluctuations  is  to  set  up  a  variation 

around  the  ideal  voltage  shape  and  since  the  ion  velocity  is  v^  ~  V*5, 
one  finds  that 


5vi  _  .  6V 
vi  H  V  ’ 


The  spread  in  arrival  times  to  the  target  is  derived  from  L  =  v.-t  where 

•  a 

5(vita)  *  0.  Thus 


it,  i! 1,  ,6V 


Here  L  is  the  target  distance  and  t.  is  the  arrival  time  on  target.  The 

a 

arrival  time  is  related  to  the  pulse  time  by  a  simple  relation  to  the  L 
required  for  voltage  ramping  by  the  amount  AV  during  the  beam  pulse  time 

^■pulse' 


Using  now  L  *  v.  t  we  find  that 

*  a 

t  ,  2V 

*a  AV  >ilse 

The  result  for  the  broadening  of  arrival  times  is  that 

at  «  iv  t 

01  AV  V  >ilse 

Since  one  typically  ramps  the  voltage  by  AV  <.  %V  and  since  the  power  multi 
plication  is  given  by  M  =  ■  one  finds  that 

1. 

The  spread  in  arrival  time  due  to  betatron  orbits  is  found  by  comparing 
motion  of  ions  on  axis  moving  at  full  vz  *  v.-  with  ions  with  finite  a 
moving  at  vz  »  (v^  -  vp*5,  where  vp  *  v^a  cos  u>gt.  For  small  a 

vz  *  v.(l  -  \  a2  cos2  Ugt)  . 

Averaging  over  one  betatron  orbit,  <  vz  >  a  v^(l  -  k  a2).  This  gives  a 
change  in  arrival  time  of  6t  *  k  a2  tg  >.  a2  tpu-|se  which  for  a  *  0.1  is 
only  0.01  of  the  original  pulse  length  adding  10%  to  a  beam  bunched  by 
M  =  10.  Even  for  a  a  0.2  this  adds  only  40%  pulse  broadening  (for  M  *  10). 
We  thus  conclude  that  power  multiplication  by  factors  of  5  is  an  easy  task 
since  only  about  10%  accuracy  in  voltage  shaping  is  required  and  no  signi¬ 
ficant  effect  is  expected  from  betatron  orbit  lengthening  for  a  <.0.2. 
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Ion  beam  propagation  in  a  fi lamented  channel  is  investigated  in 
Appendix  2.  The  fllamentation  is  a  result  of  either  radial  or  azimuthal 
current  bunching  produced  respectively  by  either  the  Weibel  or  Whistler 
instability.  In  both  cases  the  results  of  analysis  set  limits  on  the 
initial  level  of  jitter  In  the  beam  current  density.  If  the  Initial  jitter 
surpasses  these  limits,  beam  transport  will  be  affected.  When  operating 
above  the  radial  current  bunching  limit,  the  beam  will  either  be  expelled 
from  the  center  of  the  channel  or  be  pinched  into  the  center  of  the  channel 
depending  on  the  phase  of  the  initial  perturbation.  When  operating  above 
the  azimuthal  current  bunching  limit,  the  beam  is  depleted  in  time  as 
increasingly  more  ions  are  lost  to  the  dense  plasma  surrounding  the  channel. 
Typically,  these  current  bunching  limits  are  not  expected  to  be  reached  so 
that  good  transport  is  anticipated.  If  higher  current  levels  than  expected 
are  reached,  only  the  tall  of  the  beam  will  be  affected  since  the  modes 
are  nonconvecting. 

Coherent  rippling  of  the  radial  beam  envelope  is  investigated  In 
Section  III  of  Appendix  3.  This  phenomenon  Is  easily  understood  by  explor¬ 
ing  the  phase  space  picture  of  the  Ion  betatron  orbits  as  defined  by  the 
Initial  injection  conditions.  In  the  absence  of  other  mechanisms,  the 
rippling  eventually  phase  mixes  out  after  many  ripple  wavelengths  due  to 
the  weak  dependence  of  the  betatron  frequency  on  the  small  spread  in  Ion 
Injection  angles.  Other  mechanisms,  such  as  voltage  ramping,  steep  gradients 
In  the  magnetic  field  profile  (80 —  rn  for  n  >  1)  and  ion  scattering  during 
trasnport  can  considerably  reduce  this  phase  mixing  distance.  In  order  to 
minimize  variations  In  the  beam  radius,  rb(z),  at  the  beginning  of  the  trans¬ 
port  channel,  the  channel  should  be  matched  as  well  as  possible  to  the 
beam.  This  is  accomplished  by  setting  the  betatron  frequency  equal  to 


v1°rr/rs  where  vi  the  beam  velocity,  am  Is  the  maximum  injection  angle, 
and  r$  is  the  beam  spot  size. 

At  the  end  of  the  transport  channel  the  phase  mixed  density  profile  for 
the  Ideal  (B0— r)  magnetic  field  profile  shows  a  strong  peak  on  axis,  varlng 
like  r~^  as  r  -*■  0.  The  addition  of  beam  angular  momentum  reduces  the  central 
density,  however  the  beam  density  profile  will  remain  peaked  on  axis  since 
minimal  azimuthal  ion  motion  Is  observed  in  shadow-box  experiments. 

Given  the  information  obtained  from  this  work,  one  can  arrive  at  a  set 
of  optimized  plasma  and  beam  parameters  for  transport  of  a  given  current 
density  beam.  Eqs.  (11)  and  (27)  of  Appendix  3  set  upper  and  lower  limits 
on  the  required  channel  current.  Eq.  (11)  states  that  the  channel  current 
must  be  large  enough  to  confine  the  beam  and  Eq.  (27)  states  that  the  channel 
current  cannot  be  so  large  as  to  produce  a  large  degree  of  beam  rippling. 

The  electromagnetic  instability  analysis  sets  limits  on  beam  quality  (Eqs. 
(13)  and  (26)  of  Appendix  2).  The  electrostatic  instability  analysis  gives 
a  lowar  limit  on  the  plasma  density  (Eqs.  (11)  and  (31)  of  Ref.  2)  where 
the  limit  depends  on  beam  parameters  and  the  plasma  electron  temperature. 

MHD  code  work  provides  Tg(np)  for  a  given  beam.  MHD  considerations  also 
provide  a  lower  limit  on  the  plasma  density  from  plasma  expansion  considera¬ 
tion  as  discussed  earlier  in  Sec.  III-l. 

The  Ideal  target  irradiation  time  and  the  Initial  beam  pulse  duration 
determine  the  required  bunching  factor  which  in  turn  specifies  the  channel 
length  for  a  given  voltage  ramp.  The  required  bunching  factor  also  sets  an 
upper  limit  on  the  spread  in  ion  beam  injection  angles  (Eq.  (14)  of 
Appendix  3)  and  a  limit  on  the  variation  of  the  actual  voltage  ramp  from  the 

3 

ideal  ramp. 


3.  Radial  Compression  of  Ion  Beams. 

The  radial  compression  of  ion  beam  current  density  was  investigated 
for  the  possibility  of  using  either  self  magnetic  fields  or  externally 
applied  ones.  In  general,  the  gain  in  current  density  is  obtained  at  the 
expense  of  turning  a  relatively  cold  beam  into  a  hot  one.  In  principle, 
for  any  given  beam  there  will  exist  a  given  structure  that  may  bring  about 
the  highest  radial  compression  possible  Into  a  given  focus.  Since  reality 
sets  limits  on  the  possibility  of constructing  such  fields  we  have  considered 
in  the  present  studies  rather  simple  geometries  and  divided  them  into  two 
classes:  self  fields  and  external  fields.  The  case  of  steep  radial  (and/or 
axial)  gradients  shows  confinement  and  radial  compression  if  currents  beyond 
those  needed  for  uniform  current  distribution  are  driven  In  the  plasma.  The 
basic  orbits  were  studied  in  Appendices  2  and  3.  Based  on  these  results, 
recent  studies  in  that  direction  may  provide  easy  experimental  ways  of 
gaining  factors  of  4  radial  area  compression;  these  studies  are  presently 
being  pursued.  We  now  turn  to  the  weak  gradient  case. 

In  order  to  Investigate  the  value  of  radial  beam  compression.  Ion  orbits 
in  a  tapered  z-discharge  channel  were  studied  (Section  IV,  Appendix  3).  A 
multiple  scales  analysis  was  used  In  order  to  properly  treat  the  slowly 
decreasing  channel  radius.  Although  the  amplitude  of  the  Ion  betatron  orbit 
does  decrease  as  the  Ion  propagates  down  the  tapered  channel,  the  amplitude 
decreases  slower  than  the  channel  radius.  Thus,  the  effectiveness  of  radial 
compression  is  limited.  When  used  in  conjunction  with  axial  bunching,  radial 
beam  compression  can  enhance  bunching  power  multiplication  by  a  factor  of  1.3. 

The  self  magnetic  field  radial  compression  was  also  studied.  In  the 
case  of  a  fully  non-current-neutralized  Ion  beam  the  Ion  orbits  were  computed 
in  the  self-consistent  magnetic  field  assuming  zero  electric  field.  The 
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ions  are  assumed  to  enter  through  an  entrance  foil  with  small  transverse 
kinetic  energy.  As  these  orbits  converge,  the  magnetic  field  along  particle 
orbits  increases  as  r"1  as  long  as  no  particle  orbit  crossing  occurs.  This 
is  justified  for  the  case  of  a  uniform  current  distribution  (see  Appendix  4), 
and  total  ion  current  much  less  than  the  Alfven  current  (IA  *  3.1  x  107  6^ 
amperes)  and  radial  compression  ratios  of  5  giving  current  density  increases 
of  25  before  the  cold  flow  breaks  down.  In  future  studies,  attempts  may  be 
made  at  a  Vlasov  formalism  that  includes  particle  orbit  crossing  in  order 
to  see  how  much  compression  may  be  obtained.  In  addition,  realistic  time 
dependent  effects,  such  as  different  current  profiles,  voltage  variations, 
inductive  electric  fields  and  finite  rise  time  of  magnetic  fields  should  be 
studied. 

When  a  converging  sheet  ion  beam  was  treated  in  a  cylindrical  geometry 
using  magnetic  field  structures  that  have  azymuthal  symmetry  but  are  anti¬ 
symmetric  with  respect  to  a  plane  perpendicular  to  the  axis  of  symmetry.  It 
was  found  to  be  very  similar  to  the  tapered  z-discharge.  The  reason  for 
that  is  the  antisymmetric  current  flow  which,  after  converging  to  the  axis.  Is 
divided  Into  two  halves,  each  going  in  opposite  directions  through  the  two 
poles.  The  B.  field  thus  goes  as  r’1  and  drops  linearly  inside  the  plasma. 

Ion  orbits  thus  converge  geometrically  like  r”"'  and  in  addition,  are  compressed 
by  the  same  effect  as  In  the  radially  compressed  beam  in  the  tapered  z-plnch. 

Future  theoretical  work  in  this  area  of  ion  beam  focussing,  transport 
and  bunching  will  expand  on  the  work  presented  here.  In  particular,  work  is 
in  progress  studying  new  plasma  channel  schemes  that  help  to  push  upwards  the 
limits  on  deliverable  ion  currents. 
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APPENDIX  1.  Electromagnetic  Instabilities  in  a 
Focused  Ion  Beam  Propagating  Through 
A  Z-Dlscharge  Plasma. 


ELECTROMAGNETIC  INSTABILITIES  IN  A  FOCUSED 
ION  BEAM  PROPAGATING  THROUGH  A 
Z -DISCHARGE  PLASMA 

I.  INTRODUCTION 

In  an  earlier  paper,1  it  was  shown  that  focused  ion  beams  for  use  in  a  pellet  fusion  device 
can  propagate  axially  down  a  z-discharge  plasma  channel  without  generating  disruptive  micro¬ 
turbulence  due  to  electrostatic  streaming  instabilities.  The  azimuthal  magnetic  field  in  the  z- 
discharge  channel  confines  the  beam  radially  as  it  propagates.  Here  the  analysis  will  be 
extended  to  study  electromagnetic  velocity-space  instabilities.  In  particular,  the  Weibel  instabil¬ 
ity  (It  •  B  -  0.  k  ■  V,  at  0)  and  the  Whistler  instability  (k  x  |  »  0,  k  ■  V,  at  0)  are  investi¬ 
gated,  where  k  is  the  wavevector,  B  is  the  azimuthal  magnetic  field  and  V,  -  Vt  e.  is  the  axial 
streaming  velocity  of  the  beam. 

The  beam-plasma  system  consists  of  a  focused  ion  beam  (typically  a  5  MeV  proton  beam 
of  SO  ns  duration,  0.S  cm  radius,  and  a  current  of  5  x  10*<4)  propagating  down  the  axis  of  a  z- 
discharge  plasma  channel.2  The  ion  beam  is  focused  at  the  entrance  to  the  plasma  channel  (see 
Figure  1)  with  velocity  components  transverse  to  r  given  by  VJVt  at  tan  $  «  1.  A  high 
plasma  density  in  the  channel  ( n p  *  101*  cm"2)  insures  good  beam  charge  neutralization.2 
Good  beam  current  neutralization  in  the  interior  of  the  beam  also  occurs,  so  that  the  tout  mag¬ 
netic  field  is  comparable  to  that  associated  with  the  preformed  channel  established  before  beam 
injection.  The  beam  current  greatly  exceeds  that  establishing  the  channel  so  the  electron  drift 
velocity  is  approximated  by  V,  t»nbV:/nr 
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Hydrodynamic  modeling  of  the  background  pluma4  shows  that  a  uniform  channel  net* 
current  model  is  appropriate  for  the  early  times  associated  with  passage  of  the  beam  front.  This 
is  because  the  low*temperature  channel  is  established  microseconds  before  beam  injection  so 
that  complete  magnetic  diffusion  occurs.  Later  in  the  ion  pulse,  expansion  of  the  beam-heated 
high-temperature  plasma  (T  =*  25-50  eV)  reduces  the  magnetic  field  strength  in  the  interior  of 
the  channel.  The  built-up  field  in  the  expanding  cylindrical  shock  wave  is  also  enhanced  by 
significant  current  non-neutralization  in  the  cool  pluma  surrounding  the  beam-heated  channel. 
The  maximum  field  strength  just  outside  the  ion-beam  radius  can  exceed  that  established  by  the 
initial  z-discharge  current  by  a  large  factor.  Thus,  at  late  times  during  beam  passage,  the  mag¬ 
netic  field  distribution  is  closely  approximated  by  a  surface-current  model. 

In  Sec.  II,  equilibrium  models  for  such  a  beam  pluma  system  will  be  described.  In  Secs. 
Ill  and  IV,  the  Wiebel  and  Whistler  instabilities  will  be  investigated.  The  conclusions  which 
can  be  drawn  from  this  work  are  summarized  in  Sec.  V. 


II.  BEAM-PLASMA  EQUILIBRIUM 


For  mathematical  convenience,  a  slab  model  will  be  used  for  the  beam-plasma  system. 
This  is  appropriate  for  the  case  at  hand  since  ions  are  injected  into  the  channel  with  small  angu¬ 
lar  momentum  so  that  the  resulting  orbital  motion  occurs  in  a  plane.  At  early  times  in  the 
pulse,  the  net  current  (nearly  equal  to  the  channel  current)  is  uniformly  distributed  acrou  the 
channel  and  flows  in  the  z  direction.  Thus,  B  -  Byif  where 


BfXla.  |x|  <  a 

By  "  Bf,  W  >  a 


(1) 


Here,  B,  is  the  peak  value  of  the  field  and  ’a”  is  the  channel  radius.  If  the  beam  distribution 


function  fb  is  written  u 


/»(v„  v,,  v,)  -  ~  *(v,)  »(v*  +  v* 


>}  — 


ILL 

m, 


+  K) 


(2) 
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where  P:  -  myt  +  eAjc  is  the  axitl  canonical  momentum,  At  is  the  vector  potential  and  nb% 
V,  and  K  are  constants,  then 

K.  1*1  <  rb\ 

“  |0.  | Jr  |  >  rb\‘ 

where  rh  -  Vf/ta„,  K  -  V}  -  V},  a*2  -  Vzt»da  and  5fl  -  eB^ifiy.  A  smoothly  falling  density 
profile  may  be  obtained  by  replacing  the  second  delta  function  in  Eq  (2)  by  a  Maxwellian  distri¬ 
bution  function.  The  distribution  function  in  Eq.  (2)  also  states  that  all  beam  ions  cross  the  axis 
|t  the  same  angle  and  traverse  the  entire  beam  radius  during  each  betatron  oscillation.  A  more 
complicated  distribution  function  could  be  used  to  model  the  small  spread  in  angles  at  which 
the  ions  cross  the  axis,  however  little  additional  information  is  obtained  for  the  effort.  For 
mathematical  convenience,  the  form  of  the  distribution  function  given  in  Eq  (2)  will  be  used 
here.  It  is  easy  to  show  that  the  fluid  velocity  is  given  by  V  -  Kter  and  that  in  order  for  the 
beam  to  be  confined  within  the  plasma  channel,  one  must  have  rb  <  a.  Furthermore,  fb  can 
be  written  in  the  more  convenient  form 

/4(v,,  v„,  8(v,)  ®(vj  -  +  tali c2),  (4) 

where  v,  s  v4  sin$,  v,  s  Vz  +  vfl  cos  4  and  vj  -  v2  +  (v,  -  K,)2.  Here,  Vt  is  associated 
with  the  average  streaming  velocity  of  the  beam  ions  and  vb  is  associated  with  the  oscillatory 
betatron  motion  of  the  beam  ions  (vfl  <  Vfi  -  <a0rb).  The  beam-ion  orbit  equations  for  this 
field  geometry  were  solved  in  Ref.  1  and  the  results  are  summarized  in  Appendix  A. 

The  distribution  function  given  in  Eqs.  (2)  or  (4)  provides  an  appropriate  discription  of 
the  ion  beam  at  early  times  in  the  pulse.  At  late  times  By  -  0  inside  the  channel  and  the  field 
is  built-up  sharply  at  the  radial  edge  of  the  beam.  Thus  at  late  times  in  the  pulse 


/*(v,.  y$,  *) 


-~«(v,)  8(vj-  Vl),\x\  <  rb 

0,  |x|  >  rb  +  8  * 


(5) 
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where  By(rb  +  8)  is  sufficiently  large  to  coniine  the  beam  and  the  sheath  is  restricted  to  a  thin 
layer  such  that  8  «  rb.  Beam  ions  move  in  straight  line  orbits  inside  the  Held-free  channel. 
Within  the  layer  y,  <  |x|  <  yb  +  8,  the  ions  reflect  off  the  magnetic  wail,  reverse  their 
transverse  velocity  and  resume  their  straight  line  trajectories  after  reentering  the  channel.  The 
distribution  function  in  Eq.  (S)  also  results  in  the  uniform  density  profile  of  Eq.  (3)  and  in  a 
fluid  velocity  given  by  V  -  V,iz  inside  the  channel. 

The  background  plasma  provides  complete  charge  and  nearly  complete  current  neutraliza¬ 
tion  of  the  beam.  In  addition,  the  plasma  also  carries  the  z-discharge  current.  The  high  density 
desired  for  good  beam  neutralization  provides  for  a  high  frequency  of  electron-plasma  ion  colli¬ 
sions,  v„  shown4  to  be  larger  than  inside  the  channel  at  all  times  during  the  pulse.  Thus,  a 
collisional  fluid  model  is  used  for  the  background  plasma  with  the  electrons  drifting  with  velo¬ 
city  Vt  =»  (nbVJnp)  et. 

III.  THE  WEIBEL  INSTABILITY 

Two  Weibel  instabilities  will  be  investigated,  the  ion  instability  and  the  electron  instabil¬ 
ity,  which  are  respectively  driven  by  the  streaming  of  the  beam  ions  and  the  electron  drift 
motion.  Lee  and  Lampe5  report  for  electron  beams  that  the  Weibel  instability  grows  at  a 
greatly  reduced  rate  when  VJ  V,  >  where  is  the  beam  plasma  frequency  and  is 

the  electron  plasma  frequency.  Molvig*  has  shown  that  is  possible  for  electron-ion  colisions  to 
restore  rapid  growth  of  the  mode.  Although  the  concern  here  is  with  ion  beams  driving  the 
instability,  again  both  beam-thermal  effects  and  collisional  effects  are  important  Ion  betatron 
motion  will  also  be  important  in  analyzing  the  ion  instability. 
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A.  Ion  Instability 


Consider  first  the  situation  late  in  the  beam  pulse,  with  fb  given  by  Eq.  (S).  For  assumed 
large  rb,  the  perturbed  distribution  function  with  E(x)  -  E  expGibc), 


fb\  “  —  J°  dr  E  4-  -^-5- 
m,  c 


(U  ^ 

— —  exp  t(kvx  -  at)  r, 

OT 


(6) 


can  be  integrated  over  r  to  give 


fb\  “ 


i4 

av. 


Arv, 


a/* 


(7) 


Here  k  -  kex  and  straight-line  unperturbed  orbits  are  used  since  B,  -  0  for  |x|  <  rb.  The 
assumption  that  njnp  «  1  allows  one  to  write  the  usual  approximate  dispersion  equation/ 
Da  =  0,  for  the  Weibel  instability  where  |u|  <  v,  and 


Da  =s  0  «  c2k 2  +  y2  +  at2y/v,  +  Ia.  (8) 

Here  at  -  iy  for  purely  growing  perturbations,6  v,  -  mtvjmn  Vz  »  Vt  (electron  drift  motion 

is  ignored  at  present)  and  lu  is  the  beam  contribution  to  Dn\ 


4 


— 4ire2 

m, 


/ 


k\2 

iy  -  kvx 


av,  ’dv. 


d\. 


(9) 


Integrating  by  parts  and  then  using  the  calculus  of  residues  to  perform  the  remaining  nontrivial 


integration  results  in 


vjbY  I  *04+2ig  2k2V'V9 

a~  k2Vl  |(1  +y2/k2Vl)U2  y 

-  kV?,V* 

~  (1  +y2/k2V})i >2 

which  reduces  to 


(10) 
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y  »  kV0 


(11) 


OTTINGER,  MOSHER,  AND  GOLDSTEIN 


for  V,  >  Vf  and  the  limits  shown.  Solving  Eq.  (8)  for  the  growth  rate,  y,  one  obtains 

0  <  k  <  k, 

y  _  |2a,V  Kfl  \x  +  2*y,ki  k  <k  <k 
y  \  »iv,  )  I1  'k’<k<k- 

2a*p*c  (2  VJ  V)  U1  (k,  -  k),  k  =*  k0 

where 

k,  -  («>,  KrV2J/I«*  Kj);  *e  -  (a»*/c)  (2  VJ  VJ>  •«. 

The  peak  growth  rate  is  given  by 


(12) 


(13) 


For  nb/np  —  10~J  and  Vp/  Vz  >  10"1,  ~  106  sec-1  at  late  times  in  the  beam  pulse  (note  that 

y'p  is  actually  overestimated  here  since  kprb  «  1).  Thus,  no  significant  growth  can  occur  since 
y'rrb  >  0.05;  here  rb  is  the  beam  pulse  length. 


At  early  times  in  the  beam  pulse,  v,  is  larger  due  to  low  channel  temperatures  so  that  one 
might  expect  the  growth  rate  to  be  larger.  It  will  be  found,  however,  that  by  including  the 
betatron  motion  of  beam  ions  in  the  analysis,  the  perturbation  is  stabilized.  In  this  case,  fb  is 
given  by  Eq.  (4),  the  ion  orbits  are  found  in  Appendix  A,  and  fbX  is  given  by 


a.  + <v-'  •  K') 


9£, 

«  bx’ 


exp(/ur) 


(14) 


It  is  now  assumed  that  EJx)  =*  Ez  cos  tor  where  k  is  restricted  to  a  discrete  set  of  values 
by  boundary  conditions  at  rb  with  krb  >  1.  This  choice  of  EJx)  is  reasonble  since  np  and  nb 
are  both  uniform  and  nb  «  n„.  Thus,  the  mode  is  expected  to  closely  resemble  the  back* 
ground  plasma  eigenmode  which  is  sinusoidal.  Here  again,  the  background  plasma  is  treated  as 
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a  collisional  fluid  since  vt  >w„.  An  approximate  algebraic  dispersion  equation  is  obtained  by 
taking  a  weighted  spatial  average  of  DaE,(x)  where  Da  is  a  differential  operator  in  x  defined  by 


+  yi  +  ^f-~ Airie  S 


dx1  v, 

and  fbX  is  given  in  Eq.  (14)  ( Ex  =  0).  The  dispersion  equation  is  then 


0  =  Da  s  -L-  f*  dx  DaE,(x)  cos  tor 
rbEt 


(15) 


Using  the  orbit  equations  for  x'  and  v'  found  in  Appendix  A  and  the  Bessel  function  iden¬ 


tities 


exp(±  iz  sin®)  -  J0(z)  +  2  £/2,(z)  cos 2 19 


/-i 


±  2/  £j2/+1(z)  sin  (2/  +  1)  9, 
i- o 


(16) 


and 


expOz  cos  9)  -  Jt{z)  +  2  J  i1  Jt(z)  cos  19.  (17) 

f-l 

the  r  integration  in  Eq.  (14)  can  be  performed.  Thus,  Eq.  (15)  takes  the  form 

0  -  c^k2  +  y2  +  ^  +  7a  (18) 

where  in  terms  of  <*i Or,  vx,  v,)  and  a2tx,  vx,  v,) 


4 


kV, 


®1  +  «2 


(19) 


The  expressions  for  a,  and  a2  are  given  in  Appendix  B.  Terms  involving  resonances  at  higher 
harmonics  of  [  i.e.  (<■»  -  m#,)"1  for  m  *  01  are  not  considered.  The  term  involving  a,  in 
Eq.  (19)  is  the  source  of  instability.  Here,  however,  because  of  the  betatron  motion,  this  term 
vanishes  and  the  mode  is  stable  at  early  times  in  the  pulse  (see  Appendix  B).  This  same  result 
is  obtained  for  odd  eigenfunctions,  Et(x)  -  Ez  sin  toe  Mathematically  the  term  vanishes 
because  the  integrand  involving  a  |  is  an  odd  function  of  x  (see  Appendix  B).  Physically  no 
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radial  current  bunching  can  occur  since  each  ion  travels  radially  back  and  forth  across  the  beam 
profile  as  it  follows  its  betatron  orbit.  However,  azimuthal  current  bunching  can  occur  as  will 
be  discussed  in  Sec.  IV. 


B.  Electron  Instability 


The  only  other  source  of  energy  available  to  drive  the  Weibel  instability  at  early  times  in 
the  pulse  is  the  drifting  electron  background.  The  drift  velocity,  however,  is  actually  subther- 
mal  since  even  before  beam  heating  occurs  V,  =»  nb  VJn ,  <«,«■(  Tfm,) tn.  Using  a  warm 
collisional  model  (v,  >  <n„  at  all  times),  the  dispersion  equation  for  the  electron- Weibel  insta¬ 
bility  is  derived  in  Appendix  C.  Setting  y  0  in  Eq.  (C7),  it  is  found  that  y  >  0  for 
0  <  k  <  k„  -  aipiVJyfi.u,ct  where  u,  -  ( T/m,)ul .  Thus  for  k  <  k0  and  >  yv„  Eq.  (C8) 
reduces  to 


o  -  i/c1+  _ btit2YHy7  +  +  *V) _ 

V,  (-y4  +  I*(y3  +  3 k2u}y2  +  2 v,k2u2y  +  2k*u,*)  ' 
where  the  ion  beam  contribution  to  Eq.  (20)  is  ignorable. 


(20) 


At  early  times  in  the  pulse,  u^/cv,  <  1  and 


■Ji  VjCV, 


0  <  k  <  k\ 
*i  <  k  <  k0 

(*,  -  k),  k  ■=  ka 


(21) 


where  -  («*/c)  («^K,/ci»,) ,/l.  Here  T  »  4  eV,  np  s*  2  x  1011  cm-3  and  njn ,  »  SxlO”4, 
so  that  the  peak  growth  rate  is  on  the  order  of  yf  -  •  4.2 r»  *. 

At  later  times  in  the  pulse,  the  beam  has  heated  the  plasma;  then  mjtjcv,  >  1  and 


(*2KfV)>'3.  0  <  *  <  *, 

V}vJ2u ?•  k\  <  k  <  kj 


(22) 
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where  /c,  -  v,  k^-v,  VJu}  and  V,  <  uh  Now  >/-  w^K,j/(23/jcm,)  =»  l-Or**  for 

T  a»  3S  tV,  11,  -  1  X  10'*  cm”3  and  **/«,  —  1  x  10-3.  Thus,  y/r*  decreases  by  a  factor  of 
four  as  the  beam  passes  through  and  heats  the  plasma.  In  fact,  m^ujcv,  —  1  at  /  =  10.0  ns 
into  the  pulse  so  that  a  total  of  about  1.6  e-folds  occur  during  the  transit  of  the  beam. 


Note  that  rb  is  the  appropriate  time  scale  for  instability  growth  since  V,  —  0  except  dur¬ 
ing  the  passage  of  the  beam.  Furthermore,  for  small  kt  (i.e.  k}  <  IwjJZc*),  it  can  be  shown 
from  Eq.  (C9)  that 


and 


3  k}V? 
2yJ 


<  V, «  vt, 


(23) 


k?V? 


(24) 


where  yt  -  Thus,  the  unstable  mode  convects  axially,  but  with  a  group  velocity,  vf, 

which  is  much  slower  than  V,  or  V,.  At  any  given  point  in  the  plasma  the  mode  grows  only  for 
a  time  of  order  rb.  For  larger  k,  and  fixed  kx  the  mode  transforms  into  the  electrosutic 
streaming  insubility  which  was  found  to  be  stable  in  Ref.  1. 


C.  Summary  for  the  Weibel  Instabilities 


In  summary,  it  is  found  that  the  betatron  motion  of  the  beam  ions  stabilizes  the  ion- 
Weibel  instability  at  the  beam  front  while  growth  is  too  slow  at  the  tail  of  the  beam  to  allow  for 
even  one  e-fold  (y^rb  >  0.05).  The  electron- Weiber  insubility,  on  the  otherhand,  grows 
fastest  (y/r4  4.2)  at  the  front  of  the  beam  where  the  plasma  is  relatively  cold.  At  the  tail  of 

the  beam  y/r4  =■  1.0.  Although  the  electron-Weibel  insubility  grows  faster  than  the  ion- 
Weibel  insubility,  it  also  is  not  expected  to  grow  to  a  level  which  could  seriously  affect  beam 
propagation. 
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IV.  THE  WHISTLES  INSTABILITY 


The  Whistler  mode  (kxBaO,  k  •  V,  -  0)  like  the  Weibel  instability  can  be  driven 
unstable  by  particle  streaming.  However,  the  wavevector  k  -  kyi,  (i.e.  ia  in  cylindrical 
geometry)  is  perpendicular  to  the  direction  of  the  betatron  motion  of  the  beam  ions.  Hence, 
unlike  the  Weibel  instability,  the  Whistler  instability  cannot  be  stabilized  by  the  beam  ion  beta¬ 
tron  motion.  However,  a  small  spread  in  v,  (angular  momemtum)  can  reduce  the  growth  rate 
significantly.  The  electron  drift  velocity  is  ignored  when  considering  the  ion- Whistler  instability 
since  V,  »  V,. 


Consider  first  the  situation  late  in  the  beam  pulse  with  f„  given  by  Eq.  (5).  If  it  is 
assumed  for  the  moment  that  rb  is  very  large,  then  the  perturbed  distribution  function  is  given 
by 


/*!  “ 


.  F 

+  T  x  B 

8/» 

m, 

c 

8v 

/(«  -  k,Vy). 


(25) 


Using  this  expression  for  /4(,  the  perturbed  current  is  easily  calculated  and  used  to  derive  the 
dispersion  equation. 


„„ .  e,t.„  „  r. »  afc + .4  -  (26) 

Again  njn,  «  1  was  used  in  deriving  Eq.  (26)  where  w  -  iy  for  purely  growing  perturba¬ 
tions  and  V,  «  V,  (electron  drift  motion  is  ignored  at  present). 


Solving  Eq.  (26)  yields 


where 


y 


k,Vt_  0  <  ky  <  k,  I 

(kWofr./oy13,  k\  <  ky  <  *2j 
t»p^Vt/c,  ky  >  ki 


yt  >  vfi>  k\  -  -  (atyt/c)  (t»ftyt/»ic)wl. 


(27) 
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The  peak  growth  rate  (for  ky>  *2  =»  25.0  cm-1)  is  given  by  y  VJc  =*  2  x  102  rb\  for 
a  beam  with  no  spread  in  \y. 


This  result  is  only  slightly  modified  when  finite  geometry  effects  and  the  betatron  motion 
of  the  beam  ions  and  included.  Proceeding  as  in  Appendix  B  with  k  -  kjtx  +  kyiy  and  fb 
given  in  Eq.  (4),  one  obtains 


Da  —  c2(k?  +  kf)  +  y2  +  — 1 —  + 


«**> .  mmim  G(Vt) 


(28) 


where  F(kxrJ  is  defined  in  Eq.  (B6).  The  quantity  G(kxrJ  is  defined  by 


G(kxrb)  =  ~  f  k’r>  cos  *[/„(*)  J}(Z) 

+  2  i  J2m(X)  J*(Z)  ,  dX  (29) 

*-l  \z-2ix) 

with  X  -  kjgX,  2(X)  -  (kxrJ7)  (1  -  XVk}r?)'n  and  N  -  r4  +  (sin  2kxrb)/2kx.  Note  that 
G(kx  sm  0)  -  1.  F  and  G  are  only  geometrical  factors  and  do  not  modify  the  structure  of  the 
dispersion  equation.  For  large  ky  the  peak  growth  rate  now  becomes 
yy  —  (vptVJc)  [G{kxrb)Vn,  which,  aside  from  the  geometrical  correction,  is  identical  to  y„ 
found  in  Eq.  (27). 


If  the  beam  has  a  small  spread  in  v,,  the  beam  can  be  modeled  by  the  distribution  func¬ 
tion 

/»  -  8(vJ  -  expf-v2/ V*),  (30) 

where  Vy  is  the  thermal  velocity  of  the  beam  ions.  Since  the  y  motion  is  unaffected  by  By,  the 
distribution  in  v,  will  be  the  same  at  all  points  inside  the  channel.  For  convenience  finite 
geometry  effects  and  the  effects  of  the  betatron  motion  of  the  beam  ions  are  ignored  here  since 
they  have  little  affect  on  the  Whistler  mode.  Substituting  the  distribution  function  of  Eq.  (30) 
into  Eq.  (25)  yields  an  expression  for  the  perturbed  distribution  function,  /M.  The  perturbed 
current  Jbi  can  then  be  calculated  and  the  dispersion  equation  derived. 
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Da  -  e2k2  +  y2  +  utly/v,  +  4 


In  Eq.  (31),  «  —  iy  and 


,  riv  I,  +  W**  VJB. 

*  A J-  '[  K/Or  -  k,v,) 


e*p(-v*/K/) 


(31) 


(32) 


Here,  the  v4  and  4  integrations  are  trivial  and  have  already  been  carried  out.  The  quantity  4 
is  easily  written  in  terms  of  the  usual  plasma  dispersion  function,  Z( {)*•  However,  here  (  is 
pure  imaginary  so  that 


4  -«j* 


(33) 


where  Vt  >  and  erf(x)  is  the  usual  error  function.9  For  y  >  kVy,  Eq.  (31)  reduces  to  Eq. 
(26),  but  for  y  <  kVy  the  dispersion  equation  becomes 


0  ■»  c2k 2  •+  y2  +  +  at# 


(34) 


From  this  it  is  found  that  the  peak  growth  rate,  yp  -  w#  Vjc,  for  the  case  with  Vy  -  0  is 
reduced  if  V2  >  2 vlupbcVJtl,p2.  Since  this  is  the  case  here 


kVt, 

(k2V2 


0  <  k  <  ki 

kt  <  k  <  k2 


2^kvjlL  !  _  — kj  <  k  <  *„ 


(35) 


where  k{  -  fJt*2V:,  kt  ■  k{  V 2/  Vj  and  k„  -  V?  w#  VtlcVy. 

The  peak  growth  rate  is  now 

y l,-{u#Vilc)  (2 ¥luphcVjV2*2)  <  u^VJc. 

where  yp  varies  in  time  as  T~il2{t)  through  v,.  The  number  of  e-folds,  6  during  the  beam 
transit  at  any  point  in  the  plasma  is  given  by 


(36) 
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where  it  is  demonstrtted  later  than  no  significant  wave  convection  occurs.  Here  A  7 
-  I7(r*)  -  7(0)1/  HO)  for  7  in  eFand  >;(0)  is  the  peak  growth  rate  at  the  front  of  the  beam. 
A  linear  rise  in  temperature,  T(t)/T(Q)  -  l  +  A  Tt/rk,  agrees  well  with  results  of  previous 
work.3  Thus 


2t;<0>  r,  I _ j - 1 

A  T  1  (l+4r)w  '  ’ 

For  7(tj)  40  eV,  7X0)  -  4  eV,  r*  -  50  ns,  VJVt  >  0.08  and  njn,  -  10"3,  8  <  1.0  e-folds. 


A  spread  in  V  then  reduces  the  growth  of  the  ion-Whistler  insubility  to  a  tolerable  level. 


Finally  it  can  also  be  shown  that  the  ion-Whistler  insubility  does  not  convect  with  the 
beam  when  >  0.  This  needs  to  be  verified  in  order  to  justify  using  r4  as  the  appropriate 
time  scale.  Taking  Eq.  (5)  for  fb  and  setting  k  —  kyi,  +  k,il  this  dispersion  equation  becomes 

0  «•  c3*;/  -  (w,  +  iy)1  +  ^  (w,  +  iy)  +  Ia,  (38) 

V  i 

where  u  -  w,  +  iy  and 

.  fiw  , .  +  k}(V?  +  2  K,Kflcos0  +  V\  COS2  A) 

'--1.1.  J*  <***)= 

Here  the  spread  in  v,  is  neglected  since  it  will  have  little  effect  on  the  axial  group  velocity  of 
the  perturbation.  For  mathematical  convenience  finite  geometry  effects  and  the  beutron 
motion  of  the  beam  ions  are  also  ignored.  The  resulu  of  the  present  calculation  is  actually  an 
upper  limit  on  the  group  velocity  since  the  beutron  motion  of  the  beam  ions  will  tend  to  wash 
out  any  disturbance  moving  axially  on  the  beam.  For  |«,  -  ktVt\  >  ktVfi  the  denominator  of 
the  integrand  in  Eq.  (39)  can  be  expanded  and  the  integration  is  trivial.  However,  if 
-  kjV:\  <  ktV$  the  integration  is  most  easily  done  using  the  calculus  of  residues.  For 

small  k,  and  <  V„  |«,  -  k2V,\  >  k,V$  is  appropriate  for  the  ion-Whistler  insubility.  Then 
integrating  Eq.  (39),  Eq.  (38)  can  be  solved  for  «,  and  y  yielding 
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y.leL(yi+  vl/2)  '/*  -  kt  V„  (40) 

-  y  (41) 

Clearly,  v,  -  b<oJbkz  <  Vz  for  kt  <  u>^J3c.  Thus  the  ion-Whistler  instability  does  convect  in 
the  axial  direction  for  kt  <  w^/3  c,  but  with  a  group  velocity  slower  than  the  beam  streaming 
velocity.  For  kt  >  w^/3  c  the  calculus  of  residues  can  be  used  to  evaluate  Eq.  (39),  however, 
the  mode  is  then  basically  an  electrostatic  two  stream  mode.  This  mode  has  already  been 
shown  to  be  stable  in  Ref.  1. 

In  summary  it  is  found  that,  as  expected,  the  betatron  motion  of  the  beam  ions  does  not 
affect  the  ion- Whistler  instability.  The  peak  growth  rate,  however,  can  be  reduced  to  a  toler¬ 
able  level  by  the  presence  of  a  spread  in  v,  (angular  momentum  in  cylindrical  geometry).  The 
spread  in  vr  known  to  be  present  at  injection  in  typical  experiments2  is  sufficient  to  reduce  the 
number  of  e-folds  to  less  than  1.0.  Furthermore,  the  mode  convects  axially  at  a  group  velocity 
less  than  Vz.  The  electron- Whistler  instability  was  not  considered,  since  it  will  have  properties 
similar  to  the  properties  of  the  electron-Weibel  instability  (discussed  in  Section  II.B)  for  such  a 
highly  collisional  plasma  (v,  >  «„). 

V.  CONCLUSIONS 

The  purpose  of  this  paper  was  to  study  electromagnetic  velocity-space  instabilities  gen¬ 
erated  by  a  focused  ion  beam  propagating  through  a  z -discharge  plasma.  In  particular,  the 
Weibel  instability  (k  •  B  -  0,  k  •  V,  «  0)  and  the  Whistler  instability  (k  x  B  =  0,  k  V,  = 
0)  were  investigated.  This  work  is  an  extension  of  the  work  in  Ref.  1,  where  electrostatic  insta¬ 
bilities  were  investigated. 

The  ion- Weibel  instability  (driven  by  the  streaming  energy  of  the  beam  ions)  is  found  to 
be  stabilized  by  the  betatron  motion  of  the  beam  ions  at  the  front  of  the  beam.  At  the  tail  of 


14 


NRL  MEMORANDUM  REPORT  40M  ! 

the  beam,  where  beam  ions  follow  straight  line  orbits,  the  growth  of  the  ion-Weibel  instability 
is  two  slow  to  allow  for  even  one  e-fold  (y/rt  =  0.05)  during  the  transit  of  the  beam.  The  I 

electron- Weibel  instability  (driven  by  the  drifting  plasma  electrons),  on  the  other  hand,  grows 
fastest  (y/T*  as  4.2)  at  the  front  of  the  beam  where  the  plasma  is  relatively  cold  (r  »  4  eV).  j 

At  the  tail  of  the  beam,  where  T  rises  to  about  25-50  eV,  y/r6  =e  j.o.  Although  the  electron-  ! 

$ 

Weibel  instability  grows  faster  than  the  ion-Weibel  instability,  it  also  is  not  expected  to  grow  to 
a  level  which  could  drastically  affect  beam  propagation.  Only  1.6  e-foids  will  occur  during  the 
transit  of  the  beam.  It  has  also  been  shown  for  0  <  k*  <  2w^J/3cJ,  where  k  -  M*  +  M». 
that  the  electron-Weibel  instability  does  convect  axially  but  at  a  group  velocity  much  less  than 
the  beam  velocity  (vf  -  3*XV  K,/ 2«£  <  V,  «  Vt).  Thus  the  appropriate  growth  period  for 
calculating  the  number  of  e-folds  is  just  the  beam  transit  time,  r  6,  and  growth  occurs  mostly  at 
the  tail  of  the  beam.  J 

i 

Because  the  plasma  is  highly  collisions!  (v,  >  t*„  at  all  times),  the  electron- Whistler  ins¬ 
tability  will  have  properties  similar  to  those  of  the  electron-Weibel  instability.  Thus  it  is  also 
not  expected  to  grow  to  a  level  which  could  drastically  affect  beam  propagation. 

The  ion-Whistler  instability,  as  expected,  is  not  stabilized  by  the  betatron  motion  of  the 
beam  ions.  The  peak  growth  rate,  however,  can  be  reduced  to  a  tolerable  level  by  the  presence 
of  a  spread  in  v,  (angular  momentum  in  cylindrical  geometry).  Furthermore,  for  k,  < 
the  mode  also  convects  axially  at  a  group  velocity  less  than  K,(v,  -  3 ktcVll<*l*).  Thus  rb  is 
again  the  correct  time  scale  and  only  1.0  e-fold  are  expected  to  occur  during  the  transit  of  the 
beam  in  typical  experiments. 

From  these  results  and  from  the  results  for  electrostatic  instabilities  in  Ref.  1,  it  can  be 
concluded  that  it  is  possible  to  propagate  a  focused  ion  beam,  appropriate  for  a  pellet  fusion 
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device,  through  a  z-discharge  plasma  channel  without  generating  significant  growth  of  microin* 
stabilities. 
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APPENDIX  B 


The  expression  for  ax  and  a2  found  in  Eq.  (19)  are 


N 


Y  J0(Y)  +  J2(Y) 


-xa^m 


Y  J0(x)  +  j2(x) 


aj “  n((v«  “  Vt)\\  J'(x)  Jt(Y)  +  2  £  (“1)"^*(Jr)  yi»m| 


\Ja{X)  Jt(  Y)  —  2J0(X)  J2(Y) 


+  2Ja(Y)  J2(X)  +  4  £  (-1)"  Jla(X)  Ju+xiY)  +  J2k(Y)  J'2,+x(X) 

»-i  A 


"ox** 


2Vt 


/«(*>  jx(y)  +  £  (-D-[/J#+I(Ar)  n  +  Ju-ivn  Ju+x O') 


where  X  -  toe,  T  -  k\Jt*a  and  Af  -  r4  +  sin  2krJ2k. 


(Bl) 


(B2) 


In  order  to  show  that 


7i "  ~~t£L  /-*  * 609  k*  S diy 


*KjV,  a/4 


a,  -0. 


(B3) 


r  avi 

first  write  v,  -  vf  sin  6  and  v2  -  K,  +  vfl  cos  A  and  then  perform  the  v,  and  4  integration 
yielding 


/.  - 


lx2) 


x  |v#  y,u)  y  +  /|«<Z>  y,/j(z) 


-x<*eyf/2(Z) 


Y  W>  +  y,(z)  . 


(B4) 
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where  Z  -  k\f/2m0.  Integrating  by  parts  in  vfl  to  remove  the  derivative  of  the  delta  function 
and  then  performing  the  v@  integration  results  in 


'i -  SX  ix “*  * '«<z>  '■««> 

-  +  W>  y  .f'-,/2(z)  Jiniz)  + 1  y_l/2(z)  y',/2(Z) 

+  rm{z)  ym (z)  +  y3/2(z)  A/2(Z)] 

-  -f  /’1/2<z)  jmu)  \  j9(x)  +  y2(z) ) 

C  1  \lZ-i(x) 

where  t(x)  —  ( krjl )  (1  —  x2/k2r£)l/2.  Since  the  integrand  is  an  odd  function  of  x,  /, 


(B5) 

0. 


The  integral 

h  -  /_*  <bc  cos  kx  f  d\  v,a2~J  m  o>£  F{kr„)  ,  (B6) 

is  more  complicated  and  the  associated  term  in  the  dispersion  relation  does  not  contribute  to 
instability;  hence  it  will  suffice  to  write  /2  in  terms  of  the  function  F(kr„)  as  defined  in  (B6). 


APPENDIX  C 


For  a  warm  colliaional  fluid,  the  continuity  equation  and  momentum  transport  equations 
(r,  a.  7*  and  v,  -  mtvt/mi)  are 


+  V  •  (».  v.)  -  0. 


„  f« 

ir+v*VT*" 


E  + 


T.  X  B 


(Cl) 


- V/ta  -  -  v*),  (C2) 

ma/«a 

for  a  -  i,e.  Linearizing  Eqs.  (Cl)  and  (C2)  for  perturbations  with  exp  /(tor  -  ut)  dependence 
and  solving  for  the  perturbed  current  results  in 


where 


J,+J< 


(C3) 


where  w  •  iy,  ( 1 
flj-y  +  >(  +  k2u?/y. 


(if 


Ixx "  I1  +  n,J 


/n,. 


-i*Krn,/yO|. 

la  “  »i  y’D,1  ' 

,  -  y  +  *v/y,  nf  -  yv,  +  k2u2  +  2v,kW/y  +  k*u*/y\ 


In  deriving  expressions  (C4)-(C6)  it  was  assumed  that  y  <  yf. 


(C4) 

(C5) 

(C6) 

and 


Here  the  electron  streaming  is  driving  the  instability  and  n,  =  n,  (unlike  «  1), 

thus  the  complete  dispersion  equation, 

n  n  -  D„D„  -  0,  (C7) 

must  be  used.  This  results  in  the  following  dispersion  equation 
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0  - 


2  .  y **,.  , 

’  +"n7|1  +  ~ 


Si i 

n* 


yJ  +  CJ*1  + 


y  « 


yvlSh 


yci i 


t  -town? 


ft:4 


(C8) 


In  order  to  study  the  axial  convection  of  the  electron  Weibel  instability,  the  dispersion 
equation  must  be  rederived  with  |/r,l  >  0.  For  mathematical  convenience  thermal  effects  are 
ignored  (T  -  0  in  Eq.  (C2)h  In  this  case,  for  large  kx%  the  dispersion  equation  is 


0  "  cJ*,J  —  «J  - 


lujo  I  I  2^  j  _  +  iv) 


ina>  v, 


where  «'  «■  a*  —  Vt  and  now  w  — 


v,  |  *'« 


1  - 


ikzv'*, 

W 


(C9) 
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Figure  1  —  Typical  ion  trajectory  in  the  confining  azimuthal  magnetic  Held  of  the  z -discharge  channel  after  entering  the 
channel  from  the  focussing  region  on  the  left.  Here  K0  is  the  speed  of  the  beam  ion.  a  in  the  injection  angle,  r,  is  the 
spot  size  at  injection  and  a  is  the  channel  radius. 


23 


APPENDIX  2.  Ion  Beam  Propagation 
FI lamented  Channel. 


ION  BEAM  PROPAGATION  IN  A  FILAMENTED  CHANNEL 


I.  INTRODUCTION 

The  production  of  Intense  focused  ion  beans'1'  has  led  to  the 
consideration  of  using  z- discharge  plasma  channels  to  transport  then 
several  meters  to  inertial  fusion  targets?.  Analysis  has  shown  that 
although  electrostatic  beam-plasma  streaming  modes  are  stable^,  electro¬ 
magnetic  microinstabilities  will  occur  during  transport  in  stich 
channels The  fastest  growing  mode  produces  current  bunching  of  the 
electrons  in  the  channel.  In  this  report  ion  beam  propagation  in  a 
fllamented  channel  Is  investigated  in  order  to  determine  the  effects  on 
radial  beam  containment  in  channels  and  on  the  radial  beam  density 
profile. 

In  Sec.  II,  radial  current  bunching  in  the  channel  is  considered 
and  in  Sec.  Ill,  azimuthal  current  bunching  is  considered.  Finally 
the  results  of  this  work  are  summarized  in  Sec.  IV. 

II.  RADIAL  CURRENT  BUNCHING  IN  THE  CHANNEL 

If  k  ■  k  er,  then  radial  current  bunching  occurs  in  the  channel 
and  the  net  current  density  can  be  modeled  by 

4  -[l(r)  +  i1  H(r-rb)  JQ(kr>]  «z  ,  (1) 


where 


H(r-.b) 


0  «  r  *  rb  l  _ 

r  >  rb  f 


r.  is  the  beam  radius  and  J  is  the  zero  order  Bessel  function 
o  o 

which  satisfies  the  wave  equation  in  cylindrical  coordinates.  For 
Not*:  Manuscript  submitted  January  2,  1980. 
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a  channel  of  radius  r 

c 

r  r  _ 

IQ  ■  2  n  Jc  j (r)  rdr  -  2  n  Jc  j  (r)  rdr 
o  o 

Is  the  channel  current  which  confines  the  ion  beam.  Here  k  s  a^/r^ 
was  used,  where  a^  is  the  A**1  zero  of  first  order  Bessel  function 
Jj.  The  second  term  on  the  right  hand  side  of  Eq.  (1)  is  then  the 
residual  current  density  arising  from  the  bunched  electron  current 
superimposed  on  the  unbunched  Ion  be£m  current.  Since  j (r)  Is 
established  in  an  Initially  cold  plasma  and  is  driven  by  a  capacitor 
bank  on  a  time  scale  much  longer  than  the  beam  pulse  duration, when 
the  beam  is  injected  into  the  now  highly  conducting  plasma  the  total 
net  current  resists  change.  Thus,  in  the  absence  of  bunching 

i  "ie  +ib  *  6z  ’  Where  ^  <-{  V 

The  magnetic  field  which  determines  the  Ion  motion  in  a  fi lamented 
channel  is  then  found  from 


1  !r(r  V  -f[j<r)*hJo<kr>] 


’  r  s  V 


Assuming  j (r)  •  JQ  (r/rc)8  inside  the  channel  (s  st  1), 


g  Jo  * 

B®  (s+1)  c  r®"1 


4n  j 

+  - -  J  L  (kr)  ,  r  s  rfe , 


where  r  >  r.  .  The  parameter  s  is  a  magnetic  field  shaping  factor  which 
CD 

is  equal  to  one  at  the  front  of  the  beam  and  Increases  as  one  moves 
toward  the  tail  of  the  beam.  The  beam  heated  plasma  in  the  channel  at 


the  tail  of  the  beam  expands  radially  carrying  the  channel  current  with 


it.  This  decreases  the  magnetic  field  (increases  s)  near  the  center  of 
the  beam  and  produces  a  sharply  rising  magnetic  field  at  the  edge  of  the 


2 


beam.  ^ 

The  total  energy  H,  the  axial  cononical  momentum  P  ,  and  the  canon- 

z 

leal  angular  momentum  Pg  are  all  constants  of  the  motion  for  the  ions. 
From  •  -  f  I^dr  and  Eq.  (3),  P^  is  found  to  be 


mi 
where  uu 


*  z  - 


.  »  2ner  1  /m.  c  (s+1).  When  j,/j  >  0,  it  is  clear  that  the 
cb  cJo  i  J1  Jo 


s+1 


is 


beam  ion  orbits  are  radially  confined  since  1  -  j^(kr)  2  0  and  r 

monotonically  increasing.  When  j./j  <0,  the  ion  orbits  are  also  radi- 

J.  Q 


ally  confined  if  |  j^/j 

1.4  and  r' 


l1  -  V*u>] 


is  sufficiently  small  since 
s+1 


[i  -  j0(ki)J 


is  a  monotonically  Increasing  function  of 


radius. 


Now  consider  this  result  in  more  detail.  Since  P, 


e 


axial  canonical  momentum  is  given  by 


mi  ' 


V  cosa  - 
o  o 


2gicbrc  r.\s+1 
<S+1)  lrc/  . 


the  minimum 


(5) 


u 

when  j,  ■  0.  Here  a  is  the  maximum  ion  injection  angle,  V  *  (2H/m  ) 
i.  o  ox 

and  r  is  the  beam  focal  spot  size  at  the  point  of  injection  into  the 
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channel.  The  beam  envelop  radius  is  then 

1 

stl 


f  /v - -) 

c  2wcbrc  \°  “i  / 

U  m 


where  r 


r.  when  z 

D 


V  for  an  ion  injected  with  P  *  P  . 
o  z  z 


the  usual  result  is  recovered 


Fa_  /x-jL)T 

uu  .r  V  / 

cb  c  \  i  o/ 

m  m 


w 

When  s  »  1 

(7) 
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When  j,  *  0,  the  minimum  axial  canonical  momentum  is  given  by 
.  2m. j.u  .(s+1)  r 

P1  -  P°  -  -  1  -  J  <kr  >  ,  (8) 

*  *  j  k2r  L  0  8  J 

o  c 

where  this  result  only  strictly  applies  for  |j1/J0|  <  krc<r8/rc)8/(s+l) , 
In  this  case  the  beam  envelope  radius  is  found  from 


(s+1) 


Vh 


<<1  this  reduces  to 


J  (kr“)  -  J  (kr  ) 

o  b  os 


so  that  for  j^  <  0  the  beam  radius  expands,  whereas  for  j^  >  0  the  beam 

radius  decreases.  For  larger  J  Jj/J0  j  Eq.  (9)  must  be  solved  numerically. 

When  <  0,  beam  ions  can  be  magnetically  trapped  and  prevented 

from  reaching  the  axis  if  j is  sufficiently  large.  Thus  the  density 

profile  could  be  depressed  on  axis.  For  a  given  ion  injected  into  the 

channel  at  a  radius  r  and  at  an  injection  angle  a  trapping  occurs  if 

o 


|h  >  Vpd-co.  a)  +  (r./r,.)**1' 

U„  >  [l-V1',,)]  2  Vrc(s+1)  (,+1)2 


When  s  -  1  and  a  «  1,  this  reduces  to 

hi.  [vi_  +  /fo\2l 

jJ  [Wo(kro)]  L2  "<*rc  Vo/  J 


£  1  , 


so  that  trapping  begins  when  the  bunched  current  density  exceeds  the  mag¬ 
nitude  of  the  channel  current  density  which  is  required  to  confine  the 
beam.  Furthermore,  Eq.  (9)  shows  that  the  equilibrium  beam  radius  is 


also  significantly  modified  if  |l1/J0|>  1* 

Thus  when  radial  current  bunching  occurs  in  the  channel,  good  beam 
propagation  is  still  expected  unless  |  j^/jQ J  exceeds  unity.  If 
>  1,  the  beam  will  be  expelled  from  the  center  of  the  channel  when  j^ 

<  0  or  the  beam  will  pinch  when  j^  >  0.  If  r^  becomes  larger  than  rc 
the  analysis  breaks  down  and  a  more  sophisticated  analysis  is  required, 
however,  it  has  been  established  here  that  it  is  desirable  to  maintain 
Jj/j0|  <  1.  For  j^  »  jQ  exp  (y  T^)  *  *ood  beam  propagation,  requires 
Jn  <  J  exp  (y  Tfe).  Here  y  the  growth  rate  of  the  instability, 
is  the  beam  pulse  duration  and  jQ  is  the  initial  level  of  noise  in  the 
electron  current.  Since  the  predominant  source  of  j  is  the  nonunlform- 
ltles  in  the  radial  profile  of  the  ion  beam  current  density  (remember 

+  J  ),  jitter  in  the  ion  beam  current  density  must  be  kept  below 

ia  ,  i>  «p  (_  v  rb) .  as) 

Jb  Jb 

4 

Typically  Y\  1*0  “  2.0  and  j  /j  ~  0.1  for  proposed  fusion  systems. 
Thus  the  jitter  in  the  beam  current  density  must  be  kept  below  1  -  4Z. 

4 

Since  the  fllamentatlon  instability  is  nonconvectlve  ,  these  results  only 

apply  at  the  tail  of  the  beam.  At  the  front  of  the  beam  the  limit  on 

j.  /j.  stated  in  Eq.  (13)  is  less  severe  and  can  be  found  by  replacing 
on  d 

beam  of  length  x^. 

III.  AZIMUTHAL  CURRENT  BUNCHING  IN  THE  CHANNEL 

A 

If  k  -  ktg  then  azimuthal  current  bunching  occurs  in  the  channel 
and  the  net  current  density  can  be  modeled  by 


nx/xb 


where  x  is  the  distance  from  the  front  of  the  beam  for  a 
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X  “  (j(*)  +j1H(r-rb)ain  m8)6z 


(14) 


where  m  »  3,4,5,....  will  be  considered.  Special  treatment  is  required 
for  m  -  1,2  which  is  not  considered  here  since  m«  3-7  in  typical  cases. * 


Again 


/ 


j(r,0)  rdr  -  2s/  j(r)r  dr 


*/j 


is  the  channel  current  which  confines  the  ion  beam  and  j  sin  m8  is  the 
residual  current  density  arising  from  the  bunched  electron  current  super¬ 
imposed  on  the  unbunched  ion  beam  current.  The  magnetic  field  which 


determines  the  ion  motion  in  a  fllamented  channel  is  then 
s 


(s+1) cr 


8-1 


arjjr 

c(m2-4) 


sin  m3  ,  r  s  r^, 


(15) 


4ajLi 


B  *  coa  m9  t  r  £  r 

r  c(m  -4)  b 

1  s— 1 


(16) 


where  again  j(r)  -  j  (r/r  )  was  assumed. 

o  c 

In  this  case  H  and  P  are  constants  of  the  motion  but  PQ  is  not. 


From  Eq.  (15)  Pz  is  found  to  be 


P  •  2uu. 
-JL  -  z - cb 

“i  (s+1) 


m 


a+l  2 

2uu  br  j^sin  m9 

r  (m2-4)j 
c  o 


0 


r  *  v 


(17) 


Solving  Eq.  (17)  for  z  shows  that  as  long  as 
ij 
1, 


■a  * 


2  . 

m  -4 


(s+1)' 


(18) 


where  r  r_ ,  all  ion  orbits  are  confined  within  the  channel.  For 
c  b 


larger 


V*„| 


ions  can  escape  from  the  channel  and  the  beam  density  is 
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1 


gradually  depleted. 

Even  for  small  the  beam  density  on  axis  decreases  in  time 

since  | Pq  |  increases  in  time  for  most  ions.  In  fact  Pg  is  constant  only 
for  those  few  ions  which  are  injected  at  0Q  ■  nn/m  (n  -  1,2,3,. .. ,2m) . 

|VJ0|«  1 


(Rl '  Y  -  -  R(Kr  -  R8 


,2.  /2"^co*  m0\ 


m,  id  . /P  r  ,  Sf 
i  cb  z  c 


(19) 


P  u>  .  /m.r  and  a  prime 
z  cb  1  c 


where  R  -  Kr,  T  »flt,  K 
signifies  differentiation  with  respect  to  T.  Thus,  to  zero  order  in  the 
small  parameter  | j  | 

(R29/  )  -  K^P.  /fl  -  const.  (20) 

o  o 

If  Pg  is  initially  zero  and  Rq  «  A  cos  (t+  $)  (implying  small  injection 
angles, ^  i.e.  a  «  1)  the  time  averaged  Increase  in  |Pg  |  can  be  expressed 

<lpel 


as 


r  V  a 
o  o 


-  t/T, 


(21) 


where  r  V  a  »  r(z«o)r(z»o)  and 
o  o 


mV  cos  m8 
o  o 


(22) 


Thus  in  a  timer,  |Pg|  increases  from  zero  to  rQVQa .  For  the  average  ion 
<  Icos  m3  |>  *  2/n,  r  «*  r./2,  a  «  a  /2  and 

I  O  1  O  D  O 

2 


T  A _ tf 


2mV 


(23) 
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Typically  r.  »  0.4  cm,  a  «*  0.2  rad,  m«s  5,  V  «  3.1  x  10  cm/aac  and 

DO  O 

8  —1 

ui  .  »  4  x  10  sec  so  that 
cb 


T  »2.3x  10 


-9 


J<A| 


sec. 


(24) 


Thus  the  beam  hollows  out  on  a  time  scale  of  T  even  for  small 
,  although  the  beam  is  still  confined  within  the  channel.  In 
order  to  prevent  this  hollowing  out  from  occurring  during  beam  propagation 


h'K 


one  needs  T  <  T  or 


nr  a  (m  -4)1 
D  o _ o 

2mL 


exp  (-YTb)  , 


(25) 


where  T,  **  L/V  is  the  beam  transit  time  in  a  channel  of  length  L  and 
b  o 

again  j  is  the  Initial  noise  level  in  the  electron  current.  Thus  as 
argued  at  the  end  of  Sec.  II  this  implies  that  the  jitter  in  the  ion  beam 
current  density  must  be  kept  below 
U-  ittt. 


Jbn 


Ji 


S_o_  (2  .  VcbV1 

“•  Jb  W) 


exp  (- YTfa  ) 


(26) 


,-4 


Typically  this  is  found  to  be  ~  1  -  5  x  10  for  proposed  fusion 

systems.  Since  this  low  level  of  noise  is  probably  not  achievable,  some 
hollowing  out  of  the  beam  is  likely  if  current  bunching  occurs,  however 


the  beam  will  still  be  confined  in  the  channel  for  small 


V*.  - 


This 


effect  will  be  most  pronounced  at  the  tail  of  the  beam  where  largest 

4 

growth  of  the  nonconvectlve  filamentation  instability  occurs. 

11  2  2 

j^/Jo|  ^  (m  -4)/(s+l)  the  beam  density  is  gradually  depleted 
as  increasingly  more  ions  become  unconfined.  This  depletion  will  occur 
in  addition  to  the  hollowing  out  of  the  beam  and  also  occurs  predominately 
at  the  tall  of  beam.  In  order  to  prevent  this  beam  density  depletion 


8 


one  needs 


Jbn 


(m  -4)  j 


(s+l)  ^ 


2  r  {-*V 


(27) 


For  m~  5,  s  ~  4,  j  /j,  ~  0.1  and  yr.  ~  1  -  2  at  the  tail  of  the  beam, 

O  D  D 

^bn^b  ^  1  “  4Z  is  required.  At  the  front  of  the  beam  the  condition  is 
much  less  severe  with  yr^  replaced  by  where  x  is  the  distance 

from  the  front  of  the  beam  for  a  beam  of  length  x^.  Eq.  (26)  can  be 
similarly  modified  for  the  front  portion  of  the  beam. 


IV.  CONCLUSIONS 

From  the  analysis  presented  here  it  can  be  concluded  that  good  beam 
transport  in  a  filament ed  z-discharge  channel  is  possible  as  long  as  the 
current  bunching  remains  below  certain  levels.  For  radial  current  bunch¬ 
ing  it  was  found  that  | j ^ / j q  J  should  not  exceed  unity,  which  implies  that 
exP  <-YTb>-  It  |yj„l  «c»eds  rnltj,  the  b™  will  b.  «p«lled 
from  the  center  of  the  channel  when  <  0  or  the  beam  will  pinch  when 
>  0. 

For  aximuthal  current  bunching  the  beam  hollows  out  on  a  time  scale 
T~  defined  in  Eq.  (23).  This  hollowing  out  occurs  for  all  values  of 
| j 3_/J0 1  but  «1H  not  reach  a  significant  level  if  t£  <  T*  (see  Eq.  (25). 

For  | j , / J  |  a  >  the  beam  density  is  also  gradually  depleted  as  in- 

»  1  °l  (a+1)2 

creasing ly  more  ions  become  unconfined. 

These  results  se:  limits  on  the  level  of  j^n  which  can  be  tolerated 
without  seriously  affecting  beam  transport  due  to  current  bunching  effects. 


■SbSKHB 


9 


V.  REFERENCES 


1.  S.  J.  Scephanakis,  G.  Cooperstein,  Shyke  A.  Goldstein,  D.  Mosher, 

F.  W.  Oliphant,  F.  L.  Sandel  and  F.  C.  Young,  in  IEEE  Conference 
Record-Abstracts .  1979  IEEE  International  Conference  on  Plasma 
Science.  Montreal,  Quebec.  Canada  (IEEE,  New  York,  1979)  p.79. 

2.  Shyke  A.  Goldstein,  D.  Mosher,  P.  F.  Ottinger  and  W.  W.  Hsing,  in 
Inertial  Confinement  Fusion  (Optical  Society  of  America,  San  Diego, 

CA,  1978)  p.  WB7-1. 

3.  P.  F.  Ottinger,  D.  Mosher  and  Shyke  A.  Goldstein,  Phys.  Fluids  22. 

332  (1979). 

4.  P.  F.  Ottinger,  D.  Mosher  and  Shyke  A.  Goldstein  Bull.  Am.  Phys.  Soc. 
23,  800  (1978);  manuscript  submitted  to  Phys.  Fluids. 

5.  D.  G.  Colombant,  Shyke  A.  Goldstein  and  D.  Mosher,  in  IEEE  Conference 
Record-Abstracts ,  1979  IEEE  International  Conference  on  Plasma  Science. 
Montreal.  Quebec,  Canada  (IEEE,  New  York,  1979)  p.  105. 

6.  P.  F.  Ottinger,  D.  Mosher  and  Shyke  A.  Goldstein,  in  IEEE  Conference 
Record-Abstracts .  1979  IEEE  International  Conference  on  Plasma  Science. 
Montreal.  Quebec.  Canada  (IEEE,  New  York,  1979)  p.  105. 


10 


APPENDIX  3.  Propagation  of  Intense  Ion  Beams  in 

Straight  and  Tapered  Z-Discharge  Plasma 
Channels. 


PROPAGATION  OF  INTENSE  ION  BEAMS  IN 
STRAIGHT  AND  TARRED  Z-OISCHARGE  PLASMA  CHANNELS 


P.  F.  Ottinger,  D.  Mosher,  and  Shyke  A.  Goldstein8^. 
Naval  Research  Laboratory 
Washington,  D.  C.  20375 


ABSTRACT 

A  preformed  z -discharge  plasma  channel  can  be  used  to  transport  focused 
ion  beams  appropriate  for  a  pellet-fusion  device.  During  transport,  the 
beam  can  be  compressed  axially  by  time-of-flight  bunching  when  appropriate 
ion  accelerating  voltage  waveforms  are  employed.  Single-particle  orbits  in 
such  channels  are  expressible  in  terms  of  simple  harmonic  functions  for  small 
ion  injection  angles.  In  this  work,  orbit  analysis  Is  used  to  investigate 
how  non-uniformities  or  tapering  of  the  channel  and  electric  fields  present 
in  the  channel  effect  radial  beam  confinement  and  power  multiplication  by 
bunching. 


a)  Present. address:  JAYCOR,  Alexandria,  Virginia  22304 


I..  INTRODUCTION 


In  order  to  utilize  intense  light  ion  beams  (e.g.,  2-10  MeV  protons  or 
deuterons)  in  a  pellet  fusion  device,  it  is  generally  considered  necessary 
to  propagate  the  beam  a  distance  of  2-10  m  from  the  acceleration  region  while 
employing  some  focusing  scheme  to  deliver  the  beam  to  the  target.  One  con¬ 
cept^  Involves  concentrating  the  beam  by  one  order-of-magni tude  in  radius  at 
a  point  20-50  cm  in  front  of  a  small -area  ion  diode  by  a  combination  of 

O 

geometric  and  self-magnetic-field  focusing  techniques.  At  the  focus,  the 
be«n  enters  a  z-discharge  plasma  channel  which  guides  it  to  the  vicinity  of 
tho  pellet.  The  beam  is  confined  to  the  discharge  channel  during  propagation 
by  the  azimuthal  magnetic  field  produced  by  the  current  driving  the  discharge 
and  is  delivered  on  target  with  no  additional  focusing.  Figure  1  illustrates 
trajectories  of  typical  ions  in  such  a  system. 

Since  good  propagation  requires  conditions  in  the  channel  to  provide 
nearly  complete  beam  charge  and  current  neutralization,^  the  beam  Ions  do 
not  experience  collective  effects.  Single  particle  orbits  are  determined 
primarily  from  the  equations  of  motion  for  a  beam  ion  in  the  channel 
azimuthal  magnetic  field.  When  the  ion  injection  angle  into  the  channel  is 
small,  these  betatron  orbits  can  be  expressed  in  terms  of  simple  harmonic 
functions  as  demonstrated  In  Sec.  II. 

Using  the  results  of  Sec.  II,  the  z-dependence  of  the  radial  beam 
density  profile  is  obtained  in  Sec.  Ill  by  following  the  single  particle 
orbits  of  an  ensemble  of  beam  ions  Injected  into  the  channel.  Both 
analytical  and  numerical  results  are  presented  for  rippling  of  the  beam 
envelope. 
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■T  7 i-  ViT 


The  possibility  of  using  axial  bunching  and  radial  compression  of  the 
beam  to  achieve  beam  power  multi plication4  has  led  to  a  study  of  ion  orbits 
in  a  tapehed  z-discharge  channel.  In  Sec.  IV,  the  method  of  multiple 

5 

scales  Is  used  to  find  the  dependence  of  the  beam  envelope  radius  on  the 
axial  position  in  a  tapered  channel.  The  limits  on  power-density  multi¬ 
plication  due  to  axial  velocity  dispersion  in  tapered  channels  is  also 
determined. 

The  effects  on  the  ion  orbits  of  small  axial  and  radial  electric 
fields  induced  by  beam  passage  through  the  channel  are  investigated  in 
Secs.  V  and  VI  respectively.  These  analyses  determine  the  range  of  field 
strengths  which  permit  good  beam  propagation  without  significant  radial 
expansion  or  slowing  down  of  the  beam. 

If  the  plasma  channel  is  subject  to  magnetohydrodynamical  instabilities, 
beam  propagation  may  be  affected..  In  Sec.  VII  this  problem  is  studied  using 
a  code  which  solves  the  equations  of  motion  for  an  ion  injected  into  a 
bumpy  channel.  The  amplitude  of  the  bump  was  randomly  varied  in  order  to 
more  closely  simulate  the  nonlinear  evolution  of  sausage  instabilities. 
Analytical  results  in  the  form  of  a  Mathieu  equation  are  presented  which 
determine  the  importance  of  particle  resonances  in  a  constant-amplitude 
bumpy  channel.  The  focusing  and  stability  properties  of  ion  orbits  in 
these  channels  are  also  investigated. 

Finally  in  Sec.  VIII,  the  results  of  this  work  are  summarized  and  some 
concluding  remarks  are  made. 
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The  z-dlscharge  plasma  channel  providing  radial  confinement  for  the 
focused  ion  beam  may  be  formed  microseconds  before  beam  injection  so  that 
complete  magnetic  diffusion  occurs.  The  high  plasma  density  of  the  channel 
(ttp  **  10  cm  )  Insures  good  charge  neutralization  of  the  beam.  Good 
beam  current  neutralization  in  the  interior  of  the  beam  also  occurs 

*  6 

because  of  the  rapid  increase  in  the  conductivity  due  to  beam  heating.  ’ 

Thus,  the  total  magnetic  field  is  comparable  to  that  associated  with  the 
performed  channel  established  before  injection.  The  net  current,  which  is 
approximately  given  by  the  channel  current,  is  about  an  order  of  magnitude 
less  thu.n  the  bare  beam  current. 

Modeling  the  channel  current  by  a  uniform  current  profile,  the  azimuthal 
magnetic  field  is  given  by 


V/rC  r<rc 

I 

Vc/r'  r>rc 


(1) 


where  r  is  the  radius  of  the  current  channel.  In  a  carefully-mounted  experl- 
c 

ment,  imperfect  focusing  is  due  primarily  to  ion-orbit  deviations  produced 

by  such  time-varying,  azimuthally  symmetric  electromagnetic  fields  in  the 

diode.  These  fields  displace  the  ions  away  from  the  geometric  point  focus. 

At  the  best-focus  location,  ions  produced  on  opposite  sides  of  the  diode 

merge  into  the  same  spatial  region  creating  a  spread  in  angles  of  ions 

entering  the  transport  region.  These  diode  fields  do  not  Impart  azimuthal 

velocity  to  the  ions  so  that  the  angular  spread  in  injection  velocities 

lies  primarily  in  the  r-z  plane  rather  than  a  cone.  Since  the  experimental 
2 

techniques  used  to  focus  the  beam  are  not  expected  to  produce  .large 
azimuthal  asymnetries  in  the  diode  fields,  angular  momentum  effects  can  be 
ignored  in  the  lowest-order  ion  motion.  Modifications  of  the  results 
presented  here  due  to  azimuthal  motion  are  presently  under  investigation. 
Soma  numerical  results  are  presented  in  Secs.  Ill  and  VII. 

The  equations  of  motion  for  an  ion  confined  within  such  a  channel 
(i.e.,  r<rc)  are  just 
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F  -  ^cblr/rc. 
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(2) 


z  -  u>  .  rr/r  ,  (3) 

CD  C 

,  ,  ? 
where  u>cb  -  eBo/m^c.  Although  it  has  been  shown  that  the  general  solutions 

of  these  equations  can  be  written  in  terns  of  elliptic  functions ,  simple 

solutions  corresponding  to  the  i/i  <<  1  case  at  hand  are  derived.  These 

analytic  forms  provide  a  basis  for  calculations  to  follow. 

Normalizing  Eqs.  (2)  and  (3) 

R"  -  -  Z'R  ,  C4) 

Z"  «  R'R  (5) 

2  2 

where  R  *  Kr,  Z  -  Kz,  T  ■  (It,  K  -  u»cb/r  V  cosCi^Q  »  ^cb^o  cos3t/rc>  "  con8t  is 

the  speed  of  the  ion,  a0is  the  angle  of  injection  into  the  channel  and  a  primed 

quantity  signifies  d/dr.  Here  K  is  the  betatron  wavenumber  and  0  is  the 

betatron  frequency  of  the  beam  ion. 

If  the  ion  enters  the  channel  at  a  small  angle  to  the  axis,  then  the 

small  magnetic  field  required  to  confine  the  ions  results  in  e  **  Kr  «1,  and 

c 

a  simple  expansion  technique  can  be  used  to  solve  Eqs.  (4)  and  (5).  Details 
of  this  technique  can  be  found  in  Ref.  5.  Writing 

R  ■  eRj  +  6  Rj  +  .  .  .  »  (6) 

2  -  Zo  +  e2i  +  c2Z2  +  .  .  .  ,  C7) 

n  -  (1  +  ®ul  +  e2  u>2  +  .  .  .)*t  ,  (8) 

and  collecting  terms  of  like  powers  in  e ,  leads  to  a  set  of  equations  which 
can  solved  to  any  desired  order  in  the  small  parameter  e .  Zn  solving  these 
equations,  care  muse  be  taken  eo  rt move  secular  terms.  The  solution  of  Eqs.  (2) 
and  (3)  is  then 
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ft 

l 


.  r_V  cosa_  . 

tan*  - )*  tana0  , 

cb  o 

rrVo  s1n*°fl  k 

r  -  (r*  ♦  -£-S - •}**  . 

“cb  cos% 


and  r(0)  *  rQ,  z(0)  ■  0,  r(0)  °  VQ  sinatfand  z(0)  *  VQ  cosct^ 


The  maximum  angle  of  injection,  c^,  in  typical  experiments  and  for 
proposed  fusion  systems* is  on  the  order  of  0.1 -0.2  radians.  Since  the  ions 
should  be  confined  within  the  channel  (jr<  rc),  Eq.  (10)  requires  that 


I0(A)£(1.57  x  107  po*  VQ/c)  0  -  r|/rj) 


-1 


01) 


where  I0  Is  the  channel  current,  p  Is  the  ratio  of  the  beam  ion  mass  to 
the  proton  mass,  and  r&  is  the  beam  spot  size  at  Injection  into  the  channel 
This  same  result  can  be  obtained  directly  from  the  conservation  of  axial 
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canonical  momentum  which  states  that 

z  *  Vfl  coso0+  (eIQ/m^r*c2)  (r2  -  rj)  .  (12) 

Equation  (11)  Is  recovered  when  z  =  VQ,  aQ3  a^,  r  *  rc  and  rQ  3  r$  (also 
1  -  coso^as  ajjj/2).  For  a  5  MeV  proton  beam  with  om  3  0.2  radians,  rs  * 

0.4  cm  and  r£  ■  0.6  cm,  Eq.  (11)  states  that  IQ  >_ 113.0  kA. 

Equation  (9)  shows  that  the  z-motlon  has  both  a  streaming  component  and  an 
oscillatory  part  and  that  an  Ion  transit  time  in  the  channel  can  be  expressed 
as 

Tt  3  VVo  +  ^ 

where  z*  Is  the  channel  length.  The  0(e2)  corrections  to  Tt  are  Important 
when  axial  bunching  is  considered,4  since  different  ions  injected  Into  the 
channel  at  the  same  time  will  arrive  at  the  target  at  different  times.  The 
spread  In  arrival  times  for  a  straight  channel, 

<“a>s  *  <V4  vo>  K  *  “cb  r!/Vc>  ■  (14) 

can  be  obtained  from  time  averaging  Eq.  (12)  over  the  betatron  oscillations, 

and  sets  a  fundamental  limit  on  axial  bunching.  However,  at  the  point  of  injection 
the  beam  pulse  duration,  Tb>  is  still  typically  much  longer  than  (At^) ,  thus 

allowing  for  considerable  axial  compression  of  the  beam  by  proper  shaping  of  the 

accelerating  voltage  waveform. 
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III.  BEAM  RIPPLING 


If  the  radius  of  the  beam  envelope  varies  with  position  along  the 
axis  of  the  beam,  the  beam  is  said  to  ripple.  In  order  to  obtain  good 
radial  confinement  of  the  beam  and  prevent  dramatic  increases  in  beam  density 
at  minima  of  the  rippled  beam  radius,  beam  rippling  should  be  kept  to  a 
minimum.  This  phenomena  arises  from  coherent  oscillatory  motion  of  beam 
particles  injected  into  a  drift  region  permeated  by  a  confining  magnetic 
*  field.  Although  beam  rippling  occurs  in  many  situations,  the  coherence 
length,  z^,  can  vary  widely,  making  the  observation  of  rippling  impossible 
in  some  instances  (zffl  <  Ar>where  A^is  the  ripple  wavelength).  Coherence  is 
lost  by  phase  mising,  which  for  the  situation  considered  is  due  to  the 
dependence  of  the  betatron  frequency  on  the  injection  angle  of  the  individual 
ions.  Ions  enter  the  channel  at  small  angles  which  vary  uniformly  from  am 
to  -am,  where  is  determined  by  the  focusing  technique.  Since  the  beam  is 
focused  at  the  entrance  to  the  channel  and  the  focal  spot  is  defined  as  the 
point  where  only,  random  motions  determine  the  beam  radius,  the  distribution 
of  Injection  angles  Is  independent  of  injection  radius. 

Recall  from  Eq.  (10)  that  the  betatron  frequency  is  given  by 

“8  *  <“cb  W"  <’  '  “cb  l1/16rcVo  -  •  05) 

where  the  analysis  has  been  restricted  to  small  Injection  angles  (o<<1). 

9 

For  a  collection  of  Ions  Injected  into  the  channel  at  z  ■  0,  the  betatron 
frequency  varies  slightly  according  to  the  square  of  the  injection  angle  and 
Injection  radius.  The  coherence  length  of  the  beam  ripple  can  then  be 
determined  by  following  the  single  particle  orbits  of  a  collection  of  ions. 
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First  Eg.  (10)  is  rewritten  as 


r(e)  -  r0  cos  »^t  +  CV0%f^)  sin  u^t  ,  (16) 

where  r(t)  is  the  radial  position  of  a  given  ion  at  time  t  and  ujD  “  $u  ' V  /r  )* 

P  CO  o  c 

The  radial  velocity  is  just 

vr(t)  *  -rJUg  sinuigt  +  Voao  cos  u^t  .  .(17) 

Note  that  again  the  small  injection  angle  limit  is  taken  and  that  r(0)  a  r 

o 

and  v(0)  -  V.sin  a  °*  V  a  (also  v  (t)  w  V  a  (t)  )  .  Also  vfl(0)  «  v.  . 
x  o  o  o  o  r  o  y  go 

A  model  for  the  experimental  distribution  function  is  given  by 


Vro 


IS  Vi  <VV  V  |,  j 

•'W  W  ‘  j  „  } 

(0  .  W>rJ 


at  z*0  which  neglects  any  small  dispersion  in  v0Qand  approximates  the 
distribution  in  vy  by  a  square  distribution.  Here  H  is  the  unit  step  function 


^ 


x<x 


x>x 


a  is  the  maximum  injection  angle,  r  is  the  beam  spot  size  and  m.V^/2  equals 
in  s  x  o 

the  accelerating  voltage.  The  beam  density  profile  at  zsO  is  found  by  Integrating 
over  velocity  space  yielding 


where  n.  •  2a  V  N  is  the  beam  density  on  axis  and  r  is  the  beam  spot  size 
o  m  o  s 

or  beam  radius  at  z*0.  With  injection  into  the  channel  at  z**0,  the  beam 
can  be  propagated  forward  by  using  expressions  for  rQ  and  found  from 
Eqs.  (16)  and  (17), 

ro  »  r  cosu^e  -  (vt/^ )  sin  w^t  ,  (21] 


Voao  *  w0r  8in  +  vr  coa  V  ' 


Using  Eq.  (10)  to  relate  z  with  t, 


(23) 


v* 


V  r 
o  e 


(1  + 


3u»  .  r* 
cb 

16r  V 
c  o 


+  *dj  .  = 

16  J  * 


kz 


where  again  tq  and  aQ  have  been  written  in  terms  of  r  and  a.  The  distribution 
function  at  a  distance  z  downstream  from  the  point  of  injection  is  then 

,  *qr  v  U)3r  vr 

fb(r,z,vQ>V,a)-N|H(-|-  sin  k  z  +  ^  cos  k  z  +  am)-H(-^~  sinkz  +  ^- 

6^|rj  v0/|r  cos  k  z  -  (v^/ui^sin  k  z|j  6(V>Vo), 

for  | r  cos  kz  -  (v^)  sin  k  z  |*  rg.  Here  V=Vq  and  vQ|r  j  -  vQo[ro|are  constants 
of  the  motion. 

Equation  (23)  shows  that  the  dependence  on  a  of  the  harmonic  functions  in 

_  j 

Eq.  (24)  is  ignorable  unless  3  k  a  z/16  >  tt/4,  which  determines  the  phase 

-  W 

mixing  distance,  z  .  Here  k  s  Gu  , /r  V  ). 

m  co  c  o 

Integrating  Eq.  (24)  over  velocity  space  in  order  to  solve  for  the  beam  ' 

—  2 

density  profile  at  z  yields  n^fr.z)  for  *«*„  “  4n/3k  or,  where  i^Cr^z)  is 
given  in  Eq-y,  (Al)  -  (A3)  of  the  Appendix. 

For  z»*n>  kz  [see  Eq.  (23)3  varies  significantly  over  the  range  of 
integration  for  a  in  performing  the  velocity  space'  integration  of  Eq.  (24). 

Thus,  the  coherent  rippling  of  the  beam  is  phase  mixed  out  and  the  beam 


cos  k  z-a 


(24) 
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density  profile  becomes  independent  of  z  with  (r,z)  now  given  by  Eq.  (A4) 
of  the  Appendix.  The  besm  envelope  or  the  outermost  radial  edge  of  the 
beam  is  given  by 

(rg  |  cos(TTz/Xr)|  +  (Voaa/wp)  |  sin  (TTz/Xr)  |  ,  *«zml 

\  *  z»*ni  ‘ 

_  2  -2  2  2  k 

where  the  phase  mixed  beam  envelope  is  r,  ■  (r  +  V  a  /Sit)  . 

b  8  o  m  p 

From  Eq.  (25)  one  finds  that  the  beam  does  not  neck  down  to  a  radius 

smaller  than  r  if 
s 

V0  Ws  *  1  (26) 

Eq.  (26)  can  also  be  written  as  a  lower  limit  on  IQ  given  by  Iq  given  by 

Io(A)  <  (1.57  x  107  nc£  VQ/C)  (r^/rj).  (27) 

Since  it  is  desirable  to  avoid  such  necking  down  to  small  radii,  a  good 
quality  beam  will  satisfy  Eq.  (27)  as  well  as  Eq.  (11).  Figure  2  illustrates 
the  rippling  of  a  typical  beam  and  Fig.  3  shows  the  dependence  of  the 
coherence  length  on  the  maximum  injection  angle.  From  Eq.  (25)  one  also 
finds  that  the  ripple  wavelength  is  given  by  Xr  s  TT/k. 

Although  the  distribution  function  (Eq.18)  used  in  this  analysis  only 
approximates  the  actual  experimental  distribution,  the  ripple  wavelength  and 
coherence  length  should  not  depend  strongly  on  the  detailed  shape  of  the 
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distribution  function.  Provided  phase  mixing  does  not  depend  sensitively 
on  the  tail  of  the  a  distribution,  it  should  still  be  possible  to  obtain  a 
reasonable  estimate  of  the  coherence  length  from  Fig.  3  for  most  distribu¬ 
tion  functions.  Of  course,  this  is  onTy  true  in  the  absence  of  any  signifi¬ 
cant  scattering  of  the  beam  ions  (e.g.,  due  to  collisions  or  “bumpiness"  in 
the  channel)  or  any  other  mechanism  which  could  significantly  enhance  phase 
mixing  in  the  beam  (e.g.,  diode  voltage  ramping  for  bunching  or  nonuniform 
net  current  profiles).  The  ripple  amplitude  on  the  other  hand,  may  depend 
on  the  details  of  the  distribution  function. 

In  order  to  obtain  numerical  results  an  ensemble  of  one  hundred  particles 
was  injected  into  a  channel  at  2  ■  0.  Particles  were  injected  from  ten 
different  radii  ranging* from  0  to  rs  and  from  ten  equally  spaced  injection 
angles  ranging  from  -om  to  +  am-  Because  of  the  cylindrical  geometry  propor¬ 
tionately  more  particles  were  injected  at  larger  radii.  The  particles  were 
propagated  forward  in  the  r-z  plane  along  their  trajectories  by  integrating 
the  equations  of  motion  using  a  Runge-Kutta  routine  and  eventually  collected 
in  radial  bins  at  the  desired  distance  downstream. 

The  results  of  this  simulation  are  shown  in  Figs..  4  and  5,  which 
illustrate  the  rippling  of  a  beam  which  is  injected  into  a  channel  with  a 
distribution  function  at  2  -  0  given  by  Eq.  (18).  Figure  4  clearly 
shows  che  beam  rippling  when 


z  «z  .  For  comparison  with  analytical  results  the  beam  density  profile  at 

z  ■  o  and  z  «  X^/2  are  shown  in  the  upper  right  hand  corner.  The  numerical 

results  at  z  ■  44.5  cm  and  51.5  cm  should  be  compared  with  the  analytic  results 

at  z  ■  A^/2  and  the  numerical  results  at  z  8  41.0  cm  and  48.0  cm  should  be 

compared  with  the  analytic  results  at  z  »  o.  The  complete  analytic  form  of 

nfa  (r,z)  for  the  injected  distribution  function  given  in  Eq.  (18)  is  rather 

complicated  and  is  presented  in  the  Appendix.  The  good  agreement  between 

numerical  and  analytic  results  indicate  the  proper  functioning  of  the  code. 

When  z  >  the  coherent  rippling  of  the  beam  is  almost  completely  lost  as 

illustrated  in  Fig.  5.  Again  for  comparison  with  analytic  theory  an  approximate 

phase  mixed  density  profile  is  shown  in  the  upper  right  hand  corner  of  Fig.  5. 

This  profile  is  obtained  from  the  results  presented  in  the  Appendix.  The  distribution 

function  in  Eq.  (18)  was  chosen  since  it  closely  resembles  the  experimental  situation. 

For  the  proposed  parameters  used  in  Figures  4  and  5,  the  ripple  wavelength  is 

-  2 

A  «  7.0  cm  and  the  coherence  length  is  z  8  4rr/3ka  8  230  cm. 
r  .mm 

Thus  far  only  motion  in  the  r-z  plane  has  been  considered  since  it  has 

been  argued  that  most  of  the  velocity  transverse  to  the  channel  axis  is  in  the 

radial  direction.  In  order  to  verify  that  azimuthal  motion  does  not  drastically 

modify  the  results,  the  extreme  case  was  also  considered  where  ions  are  injected 

into  the  channel  with  equivalent  spreads  in  v  and  vQ.  Numerical  simulations 

show  results  that  are  similar  to  those  presented  in  Figs.  4  and  5,  however, 

the  addition  of  angular  momentum  somewhat  reduces  the  density  peak  on  axis. 

Beam  rippling  still  occurs  for  z  <  zm  and  is  phase  mixed  out  when  z>zn»  where 
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z  is  still  given  by  4n/3k  a  . 

m  m 


If  the  net  current  profile  is  not  uniform  so  that  the  azimuthal 

magnetic  field  does  not  rise  linearly  with  radius,  then  the  phase  mixing 

distance  can  be  considerably  reduced.  For  Bg  ~  r  ,  numerical  simulations  show 

that  coherence  is  lost  after  only  a  few  ripples  or  about  0.5  m  for  a  typical 

5  MeV  proton  beam  with  a  ■  0.2  rad.  This  enhanced  phase  mixing  is  illustrated 

n 

In  Fig.  '6.  Such  a  magnetic  field  profile  could  arise  if  beam  injection 
occurs  before  magnetic  diffusion  into  the  z-discharge  channel  is  completed. 
Enhanced  phase  mixing  might  also  occur  at  the  tail  of  the  beam  as  the  now 
highly  conducting  plasma  (and  frozen-in  field  lines)  expands  due  to  beam  heating 
of  the  channel. 

The  results  of  this  analysis  show  that  it  is  possible  for  ion  beams  to 
exhibit  coherent  rippling  over  distances  of  a  few  meters  before  phase  mixing 
occurs  provided  the  azimuthal  magnetic  field  rises  linearly  with  radius.  Since 
both  the  beam  radius  and  density  on  axis  vary  with  z,  beam  rippling  should  be 
an  important  consideration  when  choosing  such  things  as  channel  parameters 
(np»  n^,  rc>  r^,  etc.)  and  target  position.  The  channel  current,  IQ,  and 
the  channel  radius,  rc,  should  be  chosen  large  enough  so  that  the  outermost 
edge  of  the  rippled  beam  should  remain  in  the  channel  in  order  to  obtain  good 
radial  confinement  of  the  beam.  On  the  other  hand,  IQ  should  not  be  so  large 
that  the  beam  density  increases  dramatically  at  minima  of  the  rippled  beam 

8 

radius.  Such  an  occurrence  could  lead  to  local  heating  and  microturbulence. 
Thus,  the  beam  should  be  matched  with  the  channel  (Eq.  27)  as  well  as  possible. 


IV.  BEAM  ION  ORBITS  IN  A  TAPERED  CHANNEL 


Although  the  average  beam  density  on  axis  may  increase  by  more  than 
a  factor  of  two  over  the  injected  beam  density  on  axis  after  the  rippling 
of  the  beam  has  phase  mixed  out  (see  Fig*  5),  other  effects  can  cause 
the  beam  to  spread  in  radius  and  decrease  the  beam  density  on  axis  (e.g., 
see  Sec.  VII).  In  this  section  ion  orbits  in  a  tapered  channel  will  be 
determined  in  order  to  ascertain  the  limits  of  using  radial  compression  of 
the  beam  to  increase  the  power  density  delivered  at  the  target. 


If  the  channel  radius  decreases  linearly  such  that  rc(z)  s  rc(l-z/L), 

then 


2I0r 

cr£  (1  -  z/l)2 


r<  rc  (1  -  z/L), 


(28) 


where  L  is  the  taper  length,  rQ  is  the  channel  radius  at  z  3  0  and  Eq.  (28) 
only  applies  for  (L-z)2»  r2.  Here  it  is  assumed  that  the  beam  does  not 
burn  its  own  channel.  The  equations  of  motion  for  an  ion  confined  within 
this  channel  are  then 


zr 

-  , 

rcn  -  Z/L)* 

u.w  r  r 

z  *  — £5 _ 

rc  (1  -  z/L)2  , 


(29) 


(30) 


where  tuc{j  3  2  el0/m^c2rc.  Changing  to  normalized  variables,  Eqs.  (29)  and 
(30)  become 
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R"  «-Z'  R/(l  -  UVL)1 


(31) 


Z"  *  R'  R/{1  -  Z/KL)2  ,  (32) 

where  the  normalization  Is  the  same  as  In  Eqs.  (4)  and  (5).  If 
e  =Kr_«l»  then  R  and  Z  can  be  expanded  as  In  Eqs.  (6)  and  (7). 

Expanding  Eqs.  (31)  and  (32)  and  collecting  terms  of  equal  order  in 
em  yields  to  lowest  order 

Z"  *  0  ,  (33) 

Rj'  =-Z;  R1/(l  -  Z0/KL)2  ,  (34) 

where  it  is  assumed  that  (1  -  ZQ/KL)  <0(1)  tf*e  10n  15  not  t0°  close 

to  the  end  of  the  channel  at  z  *  L). 

The  solution  to  Eqs.  (33)  and  (34)  in  terms  of  the  real  variables  is 
z(t)  *  VQt  +  0(e2)  ,  (35) 

r(t)  -  r0(l  -  t/T)*5  cos  [-nun  (1  -  t/T)]  (36) 

♦  (VQ  sin  a/n)  (1  -  t/T)15  sin  [-  flTZn  (1  -  t/T)]+  0(e2)  )  * 

where  T  *  l/VQ,  fl2  «  «cb  VQ/rc,  r( 0)  »  rQ,  r(0)  =  VQ  ao>z(0)  »  VQ, 

o«l  and  (flT)“^  «1.  Using  conservation  of  energy 
0 

z  »  (VJ  -  r2)*5  *  V0  -  r2/(2  v0)  +  0(e3)  ,  (37) 

The  0(c2)  correction  to  z(t)  can  be  calculated. 
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Thus  Ion  trajectories  converge  radially  but  not  as  quickly  as  the  tapered 
channel  since  the  perpendicular  energy  also  Increases  In  time  at  the  expense 
of  the  parallel  energy.  From  Eqs.  (36)  and  (37)  one  finds  that  the  radial 
envelope  of  the  Ion  orbit  decreases  as 


f(W)  *  (i  +  '~s-  )H  O-2/L)15  , 

cb 


(38) 


and  the  angle  with  which  the  ion  crosses  the  axis  increases  as 


dr  « 

(fir  I 


a  (a2  .  rQ  ^Cb 
r(t)=0  (oo  rcVQ  cosc^ 


(*9) 


As  the  ion  approaches  z  3  L,  the  channel  acts  as  a  mirror  and  reflects  the 
ion  back  toward  z  3  0.  The  approximate  position  of  the  mirror  can  be  found 
by  setting  dr/dzlr(t)=(f  */2  in  Eq.  (39).  Because  the  analysis  breaks 
down  in  the  mirror  region,  the  orbit  equations  do  not  describe  this  bounce. 
The  analysis  is  not  extended  here  to' treat  this  region  since^for  reasons 
which  follow^a  target  would  be  placed  in  front  of  the  mirror  point. 

Since  the  current  density  scales  as 
J -{*<*)]  2  ~0  -  z/l)’1  , 

the  current  density  can  be  considerably  enhanced  by  radial  compression  of  the 
beam.  However,  if  radial  compression  is  used  in  conjunction  with  axial 
bunching,  then  the  spread  in  arrival  times  at  the  target  location  sets  a 
limit  on  the  power  multiplication  factor.  Since  radial  compression  produces 
a  spread  in  z,  optimizing  the  power  multiplication  will  in  turn  limit  the 
degree  of  radial  compression  which  should  be  used. 


where  zt  Is  the  target  position.  If  the  channel  were  not  tapered,  Eq. 

(43)  would  reduce  to  (Ata)$  given  in  Eq.  (14).  Clearly,  the  tapered 

channel  can  significantly  increase  At.,  which  Is  an  Important  consideration 

© 

when  axial  bunching  Is  employed  along  with  radial  compression. 

Under  ideal  circumstances,  the  diode  voltage  can  be  programmed  to 
bunch  the  beam  at  the  target,  the  tall  of  the  beam  catching  up  to  the  front 
of  the  beam  at  z»zt.  The  power  multiplication  factor  Is  then  expressed  as 

h  .  is.  r*vv  °±y. 

1- v  <««> 

where  Is  the  beam  pulse  length  and  r  Is  defined  In  Eq.  (38).  Here, 
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rb/Ata  measures  the  power  multiplication  due  to  axial  bunching  and  the 
expression  In  the  square  brackets  measures  the  power  multiplication  due  to 
radial  compression.  Since  this  analysis  only  applies  for  r(r$,  am,  z)  < 
r  (1-z/L)  (I.e.,  all  ions  are  confined  within  the  channel),  the  most  Interest 
ing  case  occurs  when  r(r$,  a^,  zt)  *  rc(l-zt/L).  In  this  case 

„  y;/«u0 

(l-zt/L)  in  [(1-Zt/L)-Jj  '  (45) 


Maximizing  M  with  respect  to  z^  yields 


M* 


8Vb'S 


a 


(46) 


Z£  *  l  [l-exp  (-h)] 


(47) 


and 


r*  -  rs  exp  (%)  (I-IJ/Iq)”1*  »  (48) 

where  must  exceed  Ii  and 
0  0 

I£(A)  -  2.6  x  107  ua*  VQ/c  .  (49) 

Comparing  this  result  with  the  power  multiplication  factor,  Ms,  obtained 
by  using  only  axial  bunching  in  a  straight  channel,  shows  that 

M*  -  2  Cexp(H)  -  1]  Ms-  a  1.3  Ms  .  (50) 

Thus  using  radial  compressions  In  conjunction  with  axial  bunching  can  enhance 
the  power  multiplication  factor  by  as  much  as  30X.  Radial  compression  will  be 
somewhat  less  effective  If  angular  motion  Is  present  In  the  beam. 


19 


V.  EFFECTS  OF  SMALL  AXIAL  ELECTRIC  FIELDS 

'Die  axial  electric  field  that  establishes  bean  current  neutralization 
by  driving  a  return  current  of  plasma  electrons  In  the  channel  also  slows 
down  the  ion  beam  and  causes  it  to  spread  radially.3  In  order  to  estimate 
the  extent  of  these  effects,  Eqs.  (4)  and  (5)  are  solved  with  the  addition 
of  a  smalt  constant  axial  electric  field  Ez  in  Eq.  (5).  In  that  case 

ZM  -  R'R  -  c\  ,  (51) 

where  &z  s  (ea/r  c"cb  vo  cosa J5  (E2/Bq).  Here  a  careful  analysis  is  required 
since  the  ion  kinetic  energy  now  decreases  in  time. 

Writing  £  *  e2x  and  n  *  (1  +  +  o^e1  +  .  .  .)x,  the  time  derivatives 

are  replaced  by 

3?  *  fj  *  <’  +  “lc  *  “2e’>  b  *  ■  '  '  ' 

3?  *  2el  ^-5^  ♦  (1  +  2  +  2  Uje1  +  e1)  ^  +  .  .  .  . 

with  R  and  Z  depending  on  both  $  and  n*  Expanding  Eqs.  (4)  and  (51)  and 
collecting  terms  in  powers  of  e  again  results  in  a  set  of  equations  which 
can  be  solved  to  any  desired  order  in  the  small  parameter  e.  Here,  however, 
the  axial  streaming  velocity,  V,  gradually  decreases  and  the  radial  envelope 
of  the  Ion  orbit,  r,  slowly  expands  on  the  5  time  scale.  The  elimination  of 
secular  terms  provides  the  set  of  equations  which  determine  V(5)  and  r(0, 
but  the  0(e*)  equation  of  Eq.  (4)  Is  needed  for  completeness.  The  resulting 
spread  in  radius  for  a  given  ion  becomes 

r(z)  »  (rj  +  rc  VQ  «J/wcb),s  (1  -  z/z^**5  ,  (52) 

where  the  radial  expansion  distance  Is 


The  slowing  down  of  the  streaming  velocity  results  in  a  decreasing  ion 
kinetic  energy  which  can  be  expressed  as 

'  v‘(«)  *  vl(0)  (1  -  z/zr)«  a  vl(0)  (1  -  2  z/zr),  (54) 

where  the  last  approximate  equality  applies  only  for  z  «z  . 

The  condition  for  the  neglect  of  the  effects  of  small  axial  electric 
fields  on  the  Ion  orbits  is  just  1  or 

eE^/iij  v;  «  1  (55) 

• 

for  am  «1.  Thus  there  is  a  lower  limit  set  on  the  channel  conductivity  in 
order  to  prevent  the  axial  electric  field,  E2  *  Jb/a»  required  for  current 
neutralization,  from  violating  the  condition  set  in  Eq.  (55).  Typically  one 
would  like  Ez  <  T.O  kV/cm  for  negligible  beam  energy  loss. 

In  addition,  it  can  be  shown  that  when  Ez  satisfies  Eq.  (55)  the 
spread  in  ion  arrival  times,  (Ata)z,  at  the  target  position  is  not  significantly 
Increased  by  the  presence  of  Ez<  In  fact,  in  that  case 

(Ata)z/(Ata)s  *  1  +  eEz  zt/m.V5  . 

Thus,  axial  bunching  is  unaffected  by  the  presence  of  an  axial  electric 
field  if  Eq.  (55)  is  satisfied. 
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VI.  EFFECTS  OF  SMALL  RADIAL  ELECTRIC  FIELDS 


The  effect  of  a  small  radial  electric  field  on  the  ion  beam  orbits 
can  be  understood  by  including  E«  *  E„  r/r_  in  Eq.  (4).  In  that  case 

lit 


R"  »  (-  V  +  e$r)  R  ,  (56) 

where  6r  *  (c2/wcb  rc  VQ  cosa J*  ( Er/BQ) .  The  basic  effect  is  to  weaken  the 
restoring  force  in  the  radial  equation  of  motion,  which  results  in  a  larger 
radius  beam.  Performing  the  analysis  of  Eqs.  (5)  and  (56)  using  the  expan¬ 
sions  shown  in  Eqs.  (6),  (7),  and  (8),  the  ion  orbits  can  again  be  solved  to 
any  desired  order  in  the  small  parameter  e.  From  this  analysis,  the  phase 
mixed  beam  envelope  Is  found  to  be 


rb  (1  +  Ef-  )  , 

b  ^  “cb  n  ° 


where  the  phase  mixed  beam  envelope  in  the  absence  of  is  defined  by 

?D*  [i  +  <rcVo  V“cb  cos  “m']*  '  (58) 

This  expression  is  obtained  from  r  by  setting  rQ  ■  r%  and  o0*  am. 

The  condition  from  Eq.  (57)  for  neglect  of  radial  electric  field  effects 
on  radial  expansion  is  then 


2  m^rJ  w2h 

E  «  — U-£>  ,  i. 
«  re 


r2  62 

92  x  10"4  (~ — 2 - )  volts/cm 

r.  tan2o 
c  m 


where  the  last  equality  applies  for  proton  beams  with  rfa  (cm),  r£  (cm)  and 
B0  (gauss)  measured  in  Gaussian  units.  For  a  5  MeV  proton  beam  with 


?b  a  rc  *  0*6  cm,  *  0.2  radians  and  BQ  =  40.0  kG,  Eq.  (59)  becomes 

E  «2  x  10*  kV/cm.  Typically,  this  condition  is  easily  satisfied  since 
r 

Ep  is  determined  only  by  magnetohydrodynamical  effects  and  thermal  effect.^ 

From  this  analysis,  one  can  also  show  that  axial  bunching  is  unaffected 
by  the  presence  of  small  radial  electric  fields. 

The  beam  ion  orbits  given  in  Eqs.  (9)  and  (10)  for  a  straight  channel 
or  in  Eqs.  (35)  and  (36)  for  a  tapered  channel  are  then  appropriate  If  the 
conditions  on  IQ,  E2  and  Er  set  forth, respectively,  in  Eqs.  (11),  (55),  and 
(59)  are  satisfied.  Under  these  conditions  the  beam  propagates  within  the 
channel  with  an  envelope  fj,  given  in  Eq.  (58)  for  a  straight  channel  and 
in  Eq.  (38)  for  a  tapered  channel  and  experiences  negligible  radial  spreading 
or  slowing  down  while  traveling  a  distance  zfc,  the  length  of  the  channel. 
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VII.  PROPAGATION  IN  A  BUMPY  CHANNEL 

Since  the  z-discharge  plasma  channel  is  formed  microseconds  before  beam 
injection,  the  possible  growth  of  sausage  type  instabilities  can  make  the 
channel  appear  bumpy  when  the  beam  is  injected.  It  is  of  value  to  understand 
how  this  bumpiness  will  affect  beam  propagation  over  distances  of  a  few 
meters.  Consider  the  situation  where  the  channel  radius  varies  sinusodially 
in  z  such  that  the  magnetic  field  is  given  by 


B 


6 


r  <  rc  (1  -  y  cos  tcz) 
r  >  rc  (1  -  f  cos  kz) 


(60) 


where  Bq  «  2Iq/ crc>  Iq  is  the  channel  current  and  2rr/iC  is.  the  wavelength,  of  the 

bumps.  If  A  «ljthen  the  magnetic  field  inside  the  channel  is  approximately 


Be  ~  (8o  r/rc^  ^  +  L  cos  KZ)  »  (61) 

and  the  radial  equation  of  motion  for  a  beam  ior.  which  is  confined  within 
r<  r_  (1  -  A/2)  becomes 

r  ■  -  (uc5  z  r/rc)  (1  +  A  cos  kz).  ^ 

Using  Eq.  (52)  and  the  energy  conservation  equation,  an-’  equation 
for  r(z)  can  be  obtained, 


£r 

dz* 


* 

C  0 


(1  +  A  COS  KZ)  , 


(63)  j 

i 

i 
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where  m.|V2/2  ■  const  is  equal  to  the  accelerating  voltage.  Writing  Eq.  (63) 
in  a  more  convenient  form 


d2R 

ds* 


[» +  (If)1] 


S/2 


R  (3  +  e  cos  2s) 


(64) 


where  R  *  icr/2,  s  ■  kz/2,  B  »  4?  ,  e  »  0Aiand  k  3  “cb/rcV  ™s 

q 

equation  closely  resembles  a  Mathieu  equation  and  can  be  solved  using  the 
method  of  multiple  scales.5  For  ter  /2~  0(e)  and  e«1 

C 


R  8  dlj  (Sq .  .)  +  e2R2  (s^,  s-j,  .  •  .)  +  .  •  »  > 


d  3  3.23 

Is  3  3sQ  e  ^  +  e  3s2  +  *  *  * 


B  « 0O  +  eS1  +  e202  +  •  •  •  * 


where  s^  «  eros.  Expanding  Eq.  (64)  and  collecting  terms  of  equal  order  in 
em,  one  obtains 


e’:  - L+S„Ri  *  0  » 

8SJ  0  1 

j2a 

e2:  — —  +  3-R,  «  -  2  — —  -  R,  fe,  +  cos  2s  )  , 


(65) 

(66) 


and  higher  order  equations.  Note  that  the  (dR/ds)2  term  in  Eq.  (64)  first 
appears  in  the  0{ej)  equation.  The  solution  of  these  equations  are  oscllla- 
tory  ’  except  nearS  *  n2  or  x»2k/n  where  n  is  a  positive  Integer. 

Since  typically  k  -  r~^  for  the  sausage  Instability,  n  ■  1  is  the  most 
important  case  to  consider.  Solving  Eq.  (65)  with  0O  »  1  yields 
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(67) 


R1  "  A<si>  exp  ^iso^  +  A  exp  <_is0* 


where  A  is  the  complex  conjugate  of  A.  Substituting  this  expression  for 
R-j  into  Eq.  (66)  and  removing  secular  terms  results  in  an  equation  for  A(s-j ) . 
Solving  this  equation  one  finally  obtains 

r(z)  a  (H  -  [a1  exp  (z/z)  __  a2  exp  (-  z/z)]  cos  (kz) 

-  (H  +01),s  [a^  exp  (z/z)  -  a2  exp  {-  z/z)]  sin  (kz)  ,  (68) 

where  |0^|  <H,  z  =  (2/icA)  (J«  -  Bp"15  and  a-j  and  a2  are  constants.  For 

|sp  >.%  »  the  solutions  are  oscillatory.  If  the  initial  conditions  for  a 

given  ion  are  such  that  a-j  /  0,  the  orbit  of  the  ion  will  exponentially  spread 
in  radius  with  an  e-folding  length  of 


z 

0 


4  (IcV 
A  1  wcb 


(69) 


where  zQ  applies  when 
zt  requires  zQ  >  zt  or 

a  <  i-  (Is!®)3* 

zt  wcb 


=  0. 


In  order  to  confine  such  ions  over  a  distance 


(70) 


This  requirement  is  very  difficult  to  satisfy  since  the  orbits  exponentially 

expand  in  radius.  For  typical  parameters,  such  as  those  used  previously 

In  Fig.  4,  Eq.  (70)  requires  that  the  channel  radius  not  vary  by  more 

than  0,45%!  Fortunately^the  resonant  region  is  narrow  and  ic~r“^  in  typical 
10 

systems  is  well  above  the  resonant  region. 

There  are  two  possible  solutions  to  this  resonance  problem:  (1)  prepare 
the  channel  such  that  k»  2  R  modes  will  be  selectively  excited  if  the  channel 
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becomes  unstable  (i.e.,  k  is  off  resonance  as  occurs  naturally),  or  (2) 

«p 

prepare  the  beam  such  that  the  ions  will  be  focused  if  unstable  ic  *  2k  modes 
are  excited  in  the  channel  (i.e.,  inject  ions  with  a2»  a^  for80  =  1). 


If  the  first  method  is  used  all  the  ions  will  have  stable  orbits  and 
will  remain  confined.  As  already  mentioned,  this  method  Occurs  naturally 
since  K^r”1  >  2k,  however,  the  possibility  of  obtaining  additional  focusing 

A* 

using  the  second  method  warrants  investigation.  If  the  second  method  is 
chosen  an  ion  entering  the  channel  at  rQ  must  have  an  injection  angle  of 


—I 

dz’ 


H  +3 


r 


2-0 


<i k  r, 


(71) 


in  order  to  exclude  unstable  orbits.  This  condition  is  found  by  setting 
a1  *  0  in  Eq.  (68). so  that 


r(z)  *  a2  exp  (-z/z)  cos  (kz  +  ir/4). 


*  (72) 


Here, the  ions  are  actually  strongly  focused  when  z t/z  >  1,  where  zt  is  the 
channel  length. 

In  order  to  ensure  that  an  ion  remains  at  least  on  a  stable  trajectory 
requires 

|a2/ax!  >  fT  exp  (zt/z)  ,  (?3) 

where  focusing  occurs  whenever  |a2/a-jj  exceeds  the  value  of  the  right  hand 
side  of  Eq.  (73).  From  this  result  one  finds  that  $  must  be  known  within  an 
accuracy  given  by 
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0 


exp  (-  kL  a/4) 


(74) 


» 


in  order  to  choose  the  proper  injection  conditions  for  focusing  Eq.  (71). 
With  $  *  4ti»ct,/rcV0ica  and  A«1  it  is  presently  not  possible  experimentally 
to  obtain  the  shot  to  shot  consistency  necessary  to  satisfy  Eq.  (74). 

Although  the  possibility  of  obtaining  additional  focusing  is  enticing, 
the  prospects  of  using  a  bumpy  channel  to  achieve  this  goal  are  not  very 
good.  Thus  method  (1), "which  avoids  particle  resonances  by  selecting 
k<tc/2,  appears  to  be  the  best  approach  to  prevent  beam  expansion.  Using 
x-r^1  and  Eq.  (11)  this  condition  on  ic  reduces  to  the  following  condition 

on  a 

m, 

°»i<0-5  (1  *  i/r£))4  •  <75) 

Thus,  for  a  given  channel  radius  and  beam  spot  size  smaller  injection  angles 
are  preferable  in  order  to  minimize  this  resonant  effect.  In  the  numerical 
results  which  follow,  the  worst  case  is  considered  with  a  -  0.2  rad  which  is 
at  the  upper  end  of  the  proposed  operating  range  of  0.1  -  0.2  rad. 

To  this  point,  the  analysis  has  only  dealt  with  small  variations  in  the 
channel  radius.  When  larger  variations  in  the  channel  radius  occur  , 
analysis  of  Mathieu's  equation  shows  that  the  unstable  resonant  region  about 
S  *  n2  broadens.9  Since  A«1  no  longer  holds,  the  previous  analysis  of 
Eq.  (63)  breaks  down  and  a  numerical  analysis  is  required.  The  equations 
of  motion  were  solved  for  an  ion  injected  into  a  bumpy  channel  where  the. 
i«roa  amplitude  bumps  were  modeled  by  a  square  waveform  which  replaces  Che 
•tnwaoidal  waveform  given  in  Eq.  (60).  Ions  move  in  Che  magnecic  field 
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2Ior/crc(z)  »  r<  pc^) 

Be  “  (76) 

2IQ/cr  ,  r>  rc(z)  , 

where  r  (z)  is  the  local  channel  radius  and  I  is  the  net  current.  The 

C  0 

amplitude  of  the  bump  was  randomly  varied  in  order  to  more  closely  simulate 
the  experimental  situation  where  the  "bumpiness"  is  not  expected  to  be  uni¬ 
form.  Ions  were  injected  into  the  channel  in  the  same  manner  as  discussed 
in  Sec.  III. 

Some  of  the  results  of  this  work  are  shown  in  Figs.  7,  8,  and  9. 

Figure  7  shows  the  percentage  of  beam  ions  which  are  lost  from  the  beam 
channel  (i.e.,  | r(zt ) | >  rc)  as  a  function  of  target  position  for  a  beam 
injected  into  a  bumpy  channel  with  a  distribution  function  at  z  =  0  given 
by  Eq.  (27).  When  2ir/<  *  3.77  cm  and  A  3  H  *  53%  of  the  beam  is  lost  at 

6.0  m  downstream  while  9%  of  the  beam  is  lost  at  6.0  m  downstream  for  A  3  h  . 

Ions  are  only  removed  from  the  system  if  |r(z)|>  2r_  at  some  point  along  the 
trajectory  of  the  ion. 

Figure  8  clearly  shows  that  there  are  fewer  ions  lost  as  the  wavelength 

of  the  bumps  is  decreased  below  about  4.0  cm.  Here)2ir/K«  6.9  cm  satisfies 

the  resonant  condition  (S  3  1  )  found  earlier  for  am  3  0.2  rad,  and  also 

m 

agrees  well  with  the  peak  of  the  curve  in  Fig.  8.  The  peak  of  the  om  3  0.1 
rad  curve  also  agrees  well  with  the  predicted  value  of  2ir/x  3  13.8  cm.  The 
resonant  wavelength  (2tr/ic)  regime  is  considerably  broadened  for  A  3  H  compared 
with  the  very  narrow  resonance  predicted  earlier  for  A«l.  The  radial 
beam  density  profile  is  plotted  as  a  function  of  axial  position  in  Fig.  9 
for  a  beam  Injected  into  a  bumpy  channel  with  Ztt/k  3  3.77  cm  and  A  3  As 
the  beam  propagates  downstream  coherence  is  quickly  lost  due  to  the  bumpiness 


29 


of  the  channel;  hencejbeam  rippling  is  unimportant  here.  As  the  beam  pro¬ 
ceeds  farther  downstream,  the  beam  gradually  expands  and  the  radial  profile 
develops  a  broad  low  density  tail.  At  6.0  m  downstream  53%  of  the  beam  ions 
are  located  at  radii  greater  than  rQt  while  20%  of  the  ions  are  located  at 
radii  less  than  rs//10  . 

Again  in  order  to  verify  that  a2imuthal  motion  does,  not  drastically  modify 

the  results,  the  extreme  case  was  considered  where  ions  are  injected  into  the 

channel  with  equivalent  spreads  in  v„  and  va.  The  numerical  results  are 

r  o 

shown  in  Fig.  10  and  should  be  compared  with  the  results  presented  in  Fig.. 

9  where  only  motion  in  the  r-z  plane  is  considered.  Since  more  particles  are 

injected  at  larger  angles  to  the  axis  when  azimuthal  motion  is  included  than 

when  only  r-z  motion  is  considered,  ions  are  lost  from  the  channel  slightly 

faster.  Approximately  75%  of  the  ions  are  lost  from  the  channel  6  m  downstream 

when  oi^  *  0.2  rad  and  2-it/k  *  3.77  cm.  As  with  case. where  only  r-z  motion 

is  considered,  much  fewer  ions  are  lost  if  om  is  reduced  to  0.1  radians  for 

m 

2k/k  *  3.77  cm. 

The  loss  rate  is  not  significantly  modified  when  nonuniform  net  current 
profiles  are  considered.  If  B0  -  r2,  numerical  simulations  show  that  60% 
of  the  ions  are  lost  from  the  channel  after  6  m  of  propagation  for  the  same 

parameters  as  used  previously  in  Fig.  9,  where  BQ  -  r  and  53%  of  the  ions 

were  lost. 

From  these  results,  it  can  be  concluded1  that,  if  possible,  modes  with 

m 

<>  2k  should  be  selectively  excited  If  the  channel  is  unstable  to  sausage 
Instability.  For  example,  if  <>.  8k  ,  then  less  than  2.0%  of  the  ions  will 
be  scattered  out  of  the  beam  in  the  first  two  meters  of  propagation  even  for 
A  as  large  as  The  most  dangerous  case  which  must  be  avoided  arises  when 


k  ■  2k  .  It  should  be  pointed  out  again  that  for  a  typical  channel ,  the  most 
unstable  sausage  mode  has  r”^  -  2.0  cm”^  which  is  off  resonance  even  for 
the  extreme  case  of  a  5  MeV  proton  beam  with  ■  0.2  rad  (2k  *  0.9  cm”*). 

For  a  5  MeV  proton  beam  with  ■  0.1  rad,  r’"*  is  well  off  resonance 
since  in  this  case  2k~  0.45  cm  . 


VIII.  CONCLUSIONS 


The  purpose  of  this  work  was  to  establish  the  groundwork  for  a  basic 
understanding  of  intense  ion  beam  propagation  in  a  z-discharge  plasma 
channel.  !t  was  implicitly  assumed  that  the  channel  remains  unaffected  by 
the  passage  of  the  beam.  Although  there  is  evidence^  that  this  may  not  be 
a  good  assumption  at  the  tail  of  the  beam,  the  analysis  presented  here  does 
correctly  describe  the  physics  of  the  forward  portion  of  the  beam.  An 
investigation  of  the  physical  processes  which  will  affect  beam  propagation 
at  the  tail  of  the  beam  is  the  subject  of  on-going  research.  Effects  on 
beam  propagation  due  to  azimuthal  beam  ion  motion  are  under  further  investigation. 


The  results  presented  here  Indicate  that  good  radial  confinement  of  the 
beam  during  propagation  In  a  straight  channel  requires  that  the  radial  and 
axial  electric  fields  remain  small.  Equations  (55)  and  (59)  set  upper  limits 
on  the  allowable  field  strengths  on  Ez  and  Er>respectively.  If  the  condition 
on  Ez  In  Eq.  (55)  is  not  satisfied,  significant  radial  expansion  and  slowing 
down  of  the  beam  is  expected  as  the  beam  propagates  downstream.  If  Er  does 
not  satisfy  Eq.  (59),  the  equilibrium  beam  radius  is  significantly  increased 
over  the  beam  radius  in  the  absence  of  Er.  These  restrictions  on  E2  and  Er 
can  be  easily  satisfied  experimentally  by  proper  choice  of  channel  parameters. 
If  the  conditions  in  Eqs.  (55)  and  (59)  are  satisfied,  it  was  also  found 
that  axial  bunching  is  unaffected  by  the  presence  of  small  radial  or  axial 
electric  fields. 

Because  of  the  small  Ion  injection  angles,  coherent  rippling  of  the 
beam  could  continue  for  distances  of  several  meters.  In  the  absence  of 
other  significant  phase  mixing  mechanisms,  the  rippling  eventually  phase 
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mixed  out  due  to  a  dependence  of  the  Ion  betatron  wavelength  on  the  Ion 
Injection  angle.  Analytic  results  show  that  the  phase-mixing  distance  Is 
given  by  zm  *  (4ir/3  a*)  (r  and  numerical  work  confirms  this  result. 

The  numerical  analysis  also  shows  that  the  beam  density  on  axis  after  phase 
mixing  can  exceed  the  beam  density  on  axis  at  injection  by  about  a  factor  of 
two.  This  corresponds  to  a  beam  radius  at  half  maximum  density  after  phase 
mixing  of  about  1//2"  times  the  original  beam  spot  size  at  Injection. 

Numerical  results  show  that  azimuthal  beam  motion  does  not  modify  the  beam  ripple 
wavelength  or  the  phase-mixing  distance  but  does  somewhat  reduce  the  final  beam 
density  on  axis. 

8eam  propagation  in  a  bumpy  channel  was  also  Investigated  in  order  to 

assess  the  problems  involved  if  the  beam  Is  injected  into  a  channel  which 

was  subject  to  growth  of  sausage  type  instabilities.  The  analysis  suggests 

that  If  the  channel  is  unstable  ic  >,  8k  modes  should  be  selectively  excited 

in  order  to  avoid  particle  resonances  which  could  deplete  the  beam  by 

radial  expansion.  For  a  typical  channel,  this  is  accomplished  with  little 

effort  since  the  fastest  growing  sausage  mode  satisfies  k  >  2k.  Using  the 

same  resonant  mechanism  to  help  focus  the  beam  does  not  seem  practical  since 
unusually  fine  control  of  the  beam  and  channel  parameters  is  required.  Numerical 

w°rh  shows  that  azimuthal  beam  motion  does  not  alter  these  results. 

Radial  compression  of  the  beam  can  be  achieved  by  injecting  the  beam 

into  a  tapered  channel.  In  this  case  the  radius  of  the  beam  converges  as 

(1  -  z/L)15  in  a  channel  which  has  a  linear  taper,  £i.e.,  r  (z)  »  r  (l-z/l)3. 

c  c 

When  used  in  conjunction  with  axial  bunching,  radial  compression  can  enhance 

the  power  multiplication  obtained  from  bunching  alone  by  as  much  as  301 
provided  the  beam  does  not  burn  its  own  channel.  Radial  compression  will  be 

somewhat  less  effective  if  angular  motion  is  present  in  the  beam. 

Future  work  will  Involve  using  the  basic  understanding  developed  here 

to  build  a  more  complete  and  detailed  picture  of  beam  propagation  in  z- 

discharge  plasma  channels. 
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An  approximate  analytic  expression  for  (r,z)  can  be  obtained  for 
the  injected  distribution  function  given  in  Eq.  (18),  if  phase  mixing 
can  be  ignored.  This  is  reasonable  as  long  as  z«zn-  In  this  case  the 
ions  can  be  propagated  forward  in  z  by  using  Eqs.  (16)  and  (17)  with 
»  nz/Xr«  Integrating  the  distribution  function  over  velocity  space 
yields 


[r2  +  r2  (z)]  /r 

0<r<r1  (z)  ] 

*)  .  2rJ(z) 

r1(z)<:r<r2(z) 

[r  +  r2  (z)]  [r3  (z)-r]  / 2r 

r2(z)<r<r3  (z) 

(X  /it)  Tan-1  (  uLr  /2V  a  ), 
r  p  a  0  m 

[r2  +  r2  (z)]  /r 

0<r<r2(z)  j 

^[r-rJ(z)]2  +  [r2(z)  +  r^z)]2^  /2r  . 

r2(z)<  r  <rx(z) 

[r  +  r2(z)]  [r3(z)-r]/2r 

r  j(z)  <  r  <r3(z) 

(X  At)  Tan"”1  (  uL  r/va)  ,  and 
r  p  s  0  tn 

|  Jtj(z)  +  r2(z)J2/r 

0  <r  <  r2(z)  | 

U  [r  -  r1(z)]2  +  [r2(z)  +  ri 2^/2r 

r2(z)  <  r<r1(z)| 

[r  +  r2  [r3*z*  ”r]  ^2r 

rj(z)<r<r3(z)  J 

for  z24z^Xx/2.  Here  r^z)  «  (V^a^/uu^ )  sin  (TTz/Xt),  r2  (z)  -  rg  cos  (n*Ar) 
-TjCxJ.rjC*)"  *8  cos  (nzAr)  +  ^  (z),  and  N  (z)  ■  2  r2/^  (z)  [r3(z)  +  ^C*)]** 

For  Xr/2^z^Xr  ,  n^  (r,z)  -  nfe  (r,  Xr-z). 

An  approximate  expression  for  the  phase  mixed  density  profile  (z  i>*m)  is 
simply  obtained  by  averaging  n^  (r,z)  over  the  ripple  length,  so  that 

.X, 

^ '  °b  Cr»*)  dz*  CM) 


This  result  is  plotted  in  the  upper  right  hand  corner  of  Fig.  5. 
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FIGURE  CAPTIONS 


1.  Trajectories  of  typical  ions  injected  into  a  z-discharge  plasma  channel 
from  the  focusing  region  to  the  left. 

2.  Illustration  of  beam  rippling  with  a  coherence  length  of  zm.  Here  the  beam 
envelope,  r^(z),  is  plotted  as  a  function  of  axial  position  for  the  distribution 

function  given  in  Eq.  (18).  For  the  proposed  parameters,  is  typically  a 

few  meters  and  r  is  about  0.5  cm. 
s 

3.  Plot  of  coherence  length'  for  beam  rippling  as  a  function  of  maximum 
injection  angle. 

4.  Illustration  of  beam  rippling  at  z<z„.  for  the  distribution  function 
given  in  Eq.  (34)  with  rc  *  0.6  cm,  wcl)  *  3.8  x  108  sec’\  VQ  =  3.1  x  109 
cm/sec,  rs  *  0.4  cm,  and  am  3  0.2  rad;  nfa  is  measured  in  arbitrary  units. 

5.  Illustration  of  nearly  complete  loss  of  coherent  beam  rippling  at 
z>  zm  3  230.6  cm  for  the  same  beam  as  shown  in  Fig,  4. 

6.  Illustration  of  nearly  complete  loss  of  coherent  beam  rippling  at 

z  «  50  cm  due  to  enhanced  phase  mixing  caused  by  a  nonuni form  channel 
current  profile.  Here  B0  ~  r*  and  all  beam  parameters  are  the  same  as 
In  Fig,  (4. 

7.  Percentage  of  beam  ions  lost  vs  axial  position  for  a  beam  Injected  into 
a  bumpy  channel  with  rc  3  0.6  cm,  r$  3  0.4  cm,  VQ  3  3.1  x  109  cm/sec, 

«  0.2,  and  IQ  3  1.2  x  10s  A. 

8.  Percentage  of  beam  ions  lost  vs  wavelength  (2tt/k)  of  the  bump  for  a 

S  MeV  proton  beam  injected  Into  a  bumpy  channel.  Here  rc  3  0.6  cm 

and  r$  3  0.4  cm  for  both  the  3  0.2  rad  (IQ  3  1*2  x  10s  A)  and  the 

o  3  0.1  rad  (I  3  3.0  x  10*  A)  curves.  The  shaded  region  indicates  the 
m  o 


proposed  operating  regime  (r£  <0.6  cm)  which  could  develop  sausage  type 
instabilities. 

9.  Plot  of  the  beam  density  profile  at  various  axial  positions  for  a  beam 
injected  Into  a  bumpy  channel.  Here^  ■  2v  rc  °  3.77  cm,  A  »  4*  and 
all  other  parameters  are  the  same  as  in  Fig.  6.  n^  is  measured  in 
arbitrary  units. 

/ 

10.  Plot  of  the  beam  density  profile  at  various  axial  positions  for  a  beam 
injected  into  a  bumpy  channel.  Here^the  ions  are  injected  into  the 
channel  with  a  spread  in  Vg  equal  to  the  spread  in  vr  and  all  other 
parameters  are  the  same  as  in  Fig.  9. 
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APPENDIX  4.  Ion  Orbits  In  a  Non-Current 
Neutralizing  Environment. 


Assume  an  ion  beam  coming  in  paraxial  under  self  Bg(r,z). 

t 


->  — . 

•  > 


2*o  &  ■  ■ 

We  now  follow  ion  orbits  assuming  that  self  increases  as  the  radius 

of  a  given  ion  decreases  because  the  current  envelope  decreases.  We  also 
assume  E  =  0. 

Instead  of  writing  down  all  equations  of  motion  we  use  energy  conserva¬ 
tion  and  the  equation  for  v?(t). 


’r  *  ^  *  jf  *  Tr  ■  <-)T  (l  -(r  *))* 


(1) 


dz 


_  eB-(r)  .  r 

J vr x  B9 * 3T *  vz  * ^  •  4  c  •  vr-,>4r  +  V*  ■  *>• 


and  since 


y. 


£dr) 


v  (1  + 


M-IO-JM#*  •  jfe  •  1 


Bn  (r/R  ). 
o 


*  V  (1  +  3  ^7*  ««  (r/Ro))  *  v(i  +  F(r» 
MV 


(2) 


1 


node*  that  F  <  0  sine*  r  <  R 

o 


We  now  use  Eq.  (2)  in  Eq.  (1)  for  vf(r). 

vr  *  (-)  v  (1  -  (1  +  F(r))2)^  =  (1  -  1  -  2F  -  F2)^ 

▼r  *=  (-v)-(-F2  -  2F)^  ■  (-)v  J2  •  |F|^  (1  - 

I(R  )*  *  2  1(V  Rfl  + 

vr  =  (-)v*1.6  — ~x —  Bn  *  R^/r  (1  -  j  — Bn  —r 


(3) 


It  Is  obvious  that  we  may  integrate  ^  *  vr 
of  Interest  to  find  the  orbits  directly.  We 


dr 

dz 


v 

r 

v 

z 


(-) 


l.6i^(r 

c 


* 


Vr 


and  find  r(t).  It  is,  however, 
now  use  Eqs.  (2)  and  O) 


<l  •-!  jr*1  vr)* 

1  •  j^fcvr 


<*> 


While  Eq.  (h)  may  be  integrated  exactly  (even  if  RQ/r  becomes  big  enough  so 

that  v  <0  the  solution  gives  a  correct  back  curved)  the  only  problem  may 
z 

be  particle  orbit  crossing  and  that  may  be  cross  checked  after  solving  for 
r(z,Ro)  for  different  R^  and  see  if  crossing  occurs.  If  it  does  then  at  the 
crossing  the  whole  formalism  fails  since  the  magnetic  field  does  not  follow 

5r  ' 

In  the  following  we  restrict  ourselves  to 

j  ^  At  R0/r  «  1  (5) 

which  implies  for  V  =  1  MV  I  =0.1  MA.  that  Bn  R^/r  ~  1  to  2  is  acceptable. 
First  we  use  Eq.  (5)  la  (fc)  end  integrate  the  resultant  equation 
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dr 

dz 


I*(R  ) 

(-)  1.6  -- ®- 


V* 


!n  ^  R  /r 
o 


(6) 


change  variable  to  y  «  0n  Ro/r  =»  r  a  Rq  e 


-y 


dr  «  (»)Ro  e"y  dy 


7  R  e”y  dy 
i»  o 
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!*(»„) 
(-)  1.6  —f- 


dz 


In  order  to  arrive  to  the  sane  z (=L)  independent  of  the  R  for  a  given 

o 

T  (say  Y  r-  1  which  gives  a  radius  reduction  by  a  factor  e). 


Y  .-y  .  1.6  I*(R  ) 

!  Z  (V) 

o  y5  Ro*V4 


(7) 


Au 


we  see  that  the  solution  to  Eq.  (7)  is  a  focus  for  all  R  if  — r - =  const. 

O  K 

O 

which  is  I  *  con:R  2  achieved  only  for  a  uniform  current  distribution.  The 
o 

distance  to  the  focus  is  given  by 


R  •  V4 
o 


7  "7  a 

z*  ~f — 5 —  '  -r  s~f- 
1.6  x  rot )  o  j® 


for  Y  >  1  the  integral  does  not  grow  fasc  and  we  anyway  are  restricted  up  to 
values  of  Y  ~  2.  We  now  need  a  tabulation  of  the  integral  as  a  function  of  Y. 


,  7 

.  1 

-1.5 

2 

2.5 

f*  «?,<* 7. 

0  7*  - 

1.56 

1.72 

1.32 

1.86 

We  use  the  value  of  Y  «  1.5  (which  gives  a  radius  reduction  by  a  factor  of  4.5 
and  an  area  reduction  of  20).  The  focal  distance  is  then 
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